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113  ABSTRACT 


TECH' OTHER  Air  Force  Cambridge  Reseach  Laboratories (LC) 

L.  G.  Hauscom  Field 

Bedford ,  Massachusetts  01730 _ 

is  abstract^  analytical  investigation  of  the  dynamic  behavior  of  tethered 

balloons  was  conducted.  This  report  covers  definition  of  tethered  balloon 
systems  and  a  study  of  stability  characteristics  and  dynamic  response  to 
wind  gusts  of  tethered  balloon  systems .  ^Balloon  systems  which  are  investiga¬ 
ted  use  the  British  3J  Barrage  Balloon,  the  Vee  Balloon  and  a  Goodyear 
Aerospace  Model  No.  1649  Single-Hull  Balloon.  The  major  tether  construction 
-''is  Columbian  Rope  Company's  N0LAR0  utilizing  prestretched  polyester  filaments. 

Three  design  altitudes,  5,000,  10,000  and  20,000  feet,  are  considered. 

The  model  for  the  tethered  balloon  system  consists  of  the  streamlined 
balloon  and  a  tether  made  tip  of  three  discrete  links.  The  derivation  of  non¬ 
linear  equations  of  motion  for  this  system  in  three  dimensions  is  presented. 

The  equations  are  linearized  for  stability  analysis  and  treated  as  uncoupled 
in  che  longitudinal  and  lateral  degrees  of  freedom.  Characteristic  equations 
which  incorporate  th<"  physical,  aerodynamic,  and  mass  characteristics  of  the 
system  are  developed  and  solved  for  the  roots  which  represent  the  frequency 
and  damping  qualifies.  The  nonlinear  equations  of  motion  are  programmed  for 
solution  or.  a  digital  computer. 

An  exploratory  investigation  to  establish  the  trends  of  dynamic 
characteristics  for  various  design  parameters  is  reported.  Design  parameters 
considered  include  balloon  shape  and  altitude,  trim  angle  of  attack,  vertical 
location  of  the  bridle  confluence  point,  net  static  lift,  tail  size,  reduced 
wind  profiles,  varying  altitude  as  encountered  in  launch  and  retrieval,  pay- 
^  wind  location  above  mean  sea  level. 
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ABSTRACT 


An  analytical  investigation  of  the  dynamic  behavior  of  tethered 
balloons  was  conducted.  This  report  covers  definition  of  tethered  balloon 
systems  and  a  study  of  stability  characteristics  and  dynamic  response  to 
wind  gusts  of  tethered  balloon  systems.  Balloon  systems  which  are  investiga¬ 
ted  use  the  British  BJ  Barrage  Balloon,  the  Vee  Balloon  and  a  Goodyear 
Aerospace  Model  No.  1649  Single-Hull  Balloon.  The  major  tether  construction 
is  Columbian  Rope  Company's  NOLAFO  utilizing  prestretched  polyester  filaments. 
Three  design  altitudes,  5,000,  10,000  and  20,000  feet,  are  considered. 

The  model  for  the  tethered  balloon  system  consists  of  the  streamlined 
balloon  and  a  tether  made  up  of  three  discrete  links.  The  derivation  of  non¬ 
linear  equations  of  motion  for  this  system  in  three  dimensions  is  presented. 
The  equations  are  linearized  for  stability  analysis  and  treated  as  uncoupled 
in  the  longitudinal  and  lateral  degrees  of  freedom.  Characteristic  equations 
which  incorporate  the  physical,  aerodynamic,  and  mass  characteristics  of  the 
system  are  developed  and  solved  for  the  roots  which  represent  the  frequency 
and  damping  qualities.  The  nonlinear  equations  of  motion  are  programmed  for 
solution  on  a  digital  computer. 

An  exploratory  investigation  to  establish  the  trends  of  dynamic 
characteristics  for  various  design  parameters  is  reported.  Design  parameters 
considered  include  balloon  shape  and  altitude,  trim  angle  of  attack,  vertical 
location  of  the  bridle  confluence  point,  net  static  lift,  tail  size,  reduced 
wind  profiles,  varying  altitude  as  encountered  in  launch  and  retrieval,  pay- 
load  location,  and  wind  location  above  mean  sea  level. 
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SECTION  I 


INTRODUCTION 


The  objective  of  this  program  is  to  investigate  the  dynamic  behavior 
of  tethered  balloons  and  in  so  doing  to  establish  design  criteria  for 
tethered  balloons,  tethers  and  payloads.  The  program  is  organized  into  three 
steps : 

(1)  Definition  of  balloon  systems  for  dynamic  analysis 

(2)  A  study  of  stability  characteristics  of  tethered  ballon 
systems 

(3)  A  study  of  dynamic  response  of  tethered  balloon  systems 
to  wind  gusts 

Reference  1  presents  a  definition  of  tethered  balloon  systems  to  be 
used  for  dynamic  studies.  These  systems  are  summarized  in  this  report,. 
Reference  2  presents  a  development  of  the  study  of  stability  characteristics 
of  the  tethered  balloon  systems.  This  report  presents  further  development  of 
the  equations  of  motion  and  solutions  for.  the  dynamic  response  of  the 
tethered  balloon  system.  Recommendations  for  the  design  of  tethered  balloon 
systems  are  presented  in  this  report. 

The  system  is  defined  as  a  balloon  tethered  at  the  end  of  a  cable 
which  is  fixed  to  a  stationary  winch.  The  tether  is  represented  by  3  straight 
links,  each  of  the  same  length.  The  links  are  considered  rigid  and  connected 
to  each  other  by  frictionless  hinges.  Goodyear  Aerospace  Corporation  has 
employed  this  method  of  representation  of  the  tethered  system  for  other  studies 
such  as  that  described  in  Reference  3. 

For  design  purposes  in  the  subject  investigation  the  following  design 
parameters  were  specified  by  AFCRL. 

Payload  and  Design  Altitude 


Payload  (lb) 

Float  Altitude  (ft  above  MSL) 

P 

h 

1500 

.5,000 

1000 

10,000 

600 

20,000 

Te  Lhers 


Two  specific  tether  constructions  )re  to  be  considered,  NOLARO  by 
Columbian  Rope  Company,  and  Amgal-Mor.itoi  AA  wire  rope  by  United  States  Steel. 
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safety  factor  of  2.0  will  be  used  with  NOLARO;  1.5  for  Ainga 1-Monitor, 
ether  design  load  will  be  based  on  a  survival  wind  of  1.3  timet  the  operationa 
wind  at  balloon  altitude. 

Wind  Profile 

The  wind  profile  as  specified  is  tabulated  below: 


SECTION  II 


BALLOON  DESCRIPTION 


The  balloon  systems  evaluated  in  this  stability  investigation  have  been 
described  in  Reference  1.  The  three  balloon  types  in  these  systems  are  the 
British  BJ  Balloon,  the  Vee-Balloon*  and  the  GAC  No.  1649  Single  Hull  Balloon 
as  depicted  in  Figures  1,  2,  and  3.  The  nominal  tethered  systems  character¬ 
istics  defined  in  Reference  1  are  summarizes  in  Tables  I  and  II. 

Static  and  dynamic  aerodynamic  characteristics  for  the  balloons  have 
been  determined  from  experimental  data  where  available  and  by  analycical 
techniques  otherwise.  The  aerodynamic  characteristics  are  presented  and 
discussed  in  Reference  2  along  with  additional  mass  and  suspension  system 
geometries  which  were  not  included  in  Reference  1. 

Design  parameters  varied  in  the  stability  investigation  include  the 
trim  angle  of  attack,  vertical  location  of  the  suspension  system  confluence 
point  below  the  balloons,  the  free-static  lift,  and  the  tail  size. 

The  trim  angle  of  attack  can  be  controlled  by  locating  the  confluence 
point  of  the  suspension  system.  The  location  of  this  point  is  established 
by  the  two  coordinates  as  shown  in  Figure  1.  The  fuselage  station  from  the 
nose  and  the  waterline  below  the  centerline  of  the  balloon  defines  this  point. 
The  trim  angle  of  attack  change  as  a  function  of  bridle  apex  point  location 
has  been  calculated  and  is  given  in  Reference  2. 

The  free-static  lift  as  used  in  this  report  is  the  excess  buoyant  lift 
provided  by  the  balloon  after  the  balloon  physical  weight  and  the  weight  of 
the  tether  in  no  wind  is  supported. 

Nominal  tail  sizes  for  each  balloon  are  depicted  in  Figures  1,  2,  and 
3.  For  the  British  BJ  Balloon,  tail  size  is  increased  and  decreased  about  the 
nominal  by  changing  the  linear  dimension''  of  each  of  the  three  tails  and  main¬ 
taining  similar  proportions.  The  inters  . tion  of  the  trailing  edge  of  the 
tails  and  the  hull  is  maintained  at  the  same  point  for  all  tail  sizes.  For  the 
Vee-Balloon  the  horizontal  tail  geometry  is  maintained.  The  vertical  tails  are 
increased  in  size  in  two  steps  for  tails  below  the  hulls.  Also  investigated 
are  the  conditions  where  vertical  tails  are  i  cated  symmetrically  above  and 
below  the  hulls.  In  addition  to  establishing  aerodynamic  characteristics  for 
these  tail  configurations  as  noted  in  Reference  2,  the  increase  in  physical 
mass  and  additional  mass  was  computed  and  used  in  the  determination  of 
stability  characteristics. 

Ascent  and  descent  studies  for  a  specific  balloon  also  incorporate  the 
change  in  mass  characteristics  due  to  air  density  changes  with  altitude. 
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Table  I.  Summary  of  Balloon/Cable  Systems 
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Trim  angle  of  attack  ^  TRIM  *s  as  ^e^ne<^  by  Table  VI  of  Reference 


Table  IX.  Summary  of  Cable  Solutions 


lallocr  Type 

Altitude 

Hull 

Tether 

B.S. 

O.D. 

wt/rt 

X 

Volume 

Type 

«t) 

(ft3) 

— ^ - - 

(lb) 

(ft) 

- - 

(lb/ft) 

(ft) 

T  Length  Total 

Tathar 
Weight 

(lb)  (ft)  (lb) 


SECTION  III 


TECHNIQUES  FOR  STUDY  OF  DYNAMIC  BEHAVIOR  OF  TETHERED  BALLOONS 


I.  MATHEMATICAL  MODELS 

A  system  of  differential  equations  was  developed  (Reference  2  and 
Appendix  A  of  this  report)  that  describes  the  motion  of  the  tethered  balloon 
in  three  dimensions.  The  degrees  of  freedom  associated  with  the  motion  are 
yaw,  pitch  and  roll  of  the  balloon  abcut  its  dynamic  mass  center,  and  pitch 
and  yaw  (lateral  rotation)  of  the  tether.  There  are  a  total  of  3  +  2N  degrees 
of  freedom  where  N  is  the  number  of  links  used  to  simulate  the  tether. 

First  consider  the  longitudinal  degrees  of  freedom.  The  dependent 
variables  shown  in  Figure  4  are  0  (pitch  of  the  balloon)  and  £r  (pitch  of 
the  "r"th  link),  where  r  is  a  particular  link.  All  angles  are  shown  positive. 


Figure  4  .  Balloon  Tether  Model  in  Longitudinal  Plane 
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In  Figure  4  VbR  is  the  relative  velocity  of  the  balloon's  center  of 
gravity  with  respect  to  the  air  and  is  the  resultant  of  the  steady  wind,  the 
wind  gust,  the  balloon  translational  motion  and  the  velocity  due  to  rotation 
cf  the  balloon  about  its  center  of  mass.  The  angle  of  attack  (a)  is  the  angle 
that  the  relative  wind  forms  with  the  longitudinal  axis  of  the  balloon. 


The  lateral  degrees  of  freedom  are  displayed  in  Figure  5  which  gives 
the  front  and  top  views  of  the  tethered  balloon.  The  lateral  degrees  of 
freedom  are:  ^(yaw  of  balloon),  0  (roll  of  balloon),  and  ar  (yaw  of 

"rthlink).  All  angles  are  shown  positive. 


Figure  5.  Balloon  Tether  Model  in  Lateral  Plane 


Pertinent  geometry  of  the  tethered  balloon  and  applied  forces  are 
identified  in  Figure  6. 

In  order  to  separate  the  equations  of  motion  into  a  longitudinal 
response  and  a  lateral  response,  it  was  further  assumed  that  the  system 
was  near  equilibrium.  This  resulted  in  a  set  of  equations  describing  the 
longitudinal  motion  which  is  coupled  only  in  the  pitching  variables  of  the 
balloon  and  the  pitching  variables  of  the  tether.  However,  the  second  set 
of  equations  for  the  lateral  motion  does  not  completely  uncouple  from  the 
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Lg  -  Gross  static  lift  (Total  weight  of  displaced  air)\  |  / 

Wg  -  Total  weight  of  balloon  including  internal  gases  *TJ 
FgYA  -  Vertical  aerodynamic  forces  acting  on  balloon  \ 

fBHA“  Horizontal  aerodynamic  forces  acting  on  balloon  J — V 

Mflg  -  Aerodynamic  pitching  moment  acting  on  balloon 
CA  -  Aerodynamic  reference  center  (at  center  of  hull  volume) 

CB  -  Center  of  buoyancy  of  total  displaced  volume 

Cm  -  Apparent  mass  center  and  dynamic  center  in  pitch  (and  yaw) 

CG  -  Center  of  gr.vity  of  total  weight  of  balloon 


Figure  6.  Balloon  Geometry  and  Applied  Forces 

longitudinal  degrees  of  freedom  because  the  equilibrium  angles  in  the  longi¬ 
tudinal  plane  are  not  zero.  Therefore,  when  solving  the  lateral  degrees  of 
freedom,  it  must  be  assumed  that  the  longitudinal  variables  remain  constant 
and  equal  to  their  equilibrium  values.  In  both  the  longitudinal  and  lateral 
cases,  the  tether  is  simulated  by  three  rigid  links  The  number  of  uncoupled 
dynamic  equations  is  four  for  the  longitudinal  response  and  five  for  the 
lateral  response. 

2.  GENERAL  STABILITY  THEORY 

The  equilibrium  configuration  of  a  tethered  ballon  can  be  defined  as 
that  position  which  demands  that  the  summation  of  all  applied  moments  equals 
zero.  The  equilibrium  is  said  to  be  stable  if,  for  any  small  disturbance, 


the  system  ultimately  returns  to  its  equilibrium  conditions.  Two  types  of 
stability  are  of  interest.  In  the  first  (statically  stable),  a  small  displace¬ 
ment  of  the  system  will  create  forces  which  tend  to  return  the  system  to  its 
equilibrium  position.  The  second  (dynamically  stable)  produces  a  motion  which 
eventually  restores  equilibrium.  If  the  motion  is  periodic,  it  is  character¬ 
ized  by  a  damped  frequency  and  a  damping  ratio.  Similar  definitions  apply 
for  statically  and  dynamically  unstable  motions.  A  third  possibility  is  for 
the  system  to  be  neutrally  stable  during  which  the  motion  neither  diverges 
nor  converges. 

It  was  necessary  during  this  study  to  develop  techniques  to  investigate 
and  understand  the  stability  of  the  system.  To  this  end,  characteristic  equa¬ 
tions  are  derived.  The  general  approach  is  as  follows: 

(1)  Derive  the  nonlinear  equations  of  motion  in  three 
dimensions  for  each  degree-of-f reedom 

(2)  Assume  the  motion  is  near  equilibrium  so  that  the 
equations  can  be  linearized  and  separated  into  a 
longitudinal  motion  and  a  lateral  motion 

(3)  Laplace  transform  the  linear  equations  from  the  time 
domain  to  the  "S"  domain  assuming  that  the  initial  con¬ 
ditions  are  zero.  This  establishes  a  matrix  equation  of 
the  following  form: 

a]  {x  (S)j  =  {o}  (1) 

where  jX  (S)  }  is  the  eigenvector  and  [A]  is  a  square  matrix  whose  elements 
are  quadratics  in  S  containing  the  physical  properties  of  the  system. 

(4)  Expand  the  determinant  of  [a]  such  that  the  characteristics 
polynominal  is  obtained. 

Each  root  of  the  characteristic  equation  represents  a  term  in  the 
general  solution  of  the  form,  A^e  S^t,  where  is  the  "i"th  root  and  A^ 
is  an  amplitude,  dependent  on  the  initial  conditions  of  the  system.  Both 
real  and  complex  roots  may  appear  where  the  complex  roots  occur  in  conjugate 
pairs,  in  general  for  "n"  degrees  of  freedom,  the  characteristic  equation 
will  yield  "2n"  roots.  Each  pair  of  complex  conjugate  roots  represents  one 
oscillatory  motion,  while  each  real  root  represents  one  aperiodic  motion. 

First  consider  an  oscillatory  system.  This  motion  is  characterized 
by  two  roots  of  the  form  Si  =  X  ±  i  Y,  where  X  and  Y  are  real  numbers 

and  i  =  J~- T.  Several  important  quantities  can  be  found  from  the  root.  The 

natural  frequency  associated  with  this  motion  is  =  n/x2  +  Y2.  The 

damping  ratio  is  £  =  — gjr — •  The  damping  frequency  is  (j,  =  b)n  \fl  -  £2  =  Y. 
It  is  also  of  interest  to  Know  the  time  to  half  amplitude  for  a  stable  root  or 
tha  time  to  double  amplitude  for  an  unstable  root.  This  quantity  can  easily 
be  found  by  considering  one  oscillatory  motion.  The  general  solution  for 
free  vibration  is 
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Z  =  Ce'  ^  w  n  c  Sin  (  ,ud  t  +  0)  (2) 

where  0  is  the  phase  angle  dependent  upon  initial  conditions 

C  is  a  constant  dependent  upon  initial  conditions 

The  second  possibility  is  an  aperiodic  motion  given  by  the  expression 

Z  =  CeXt  (3) 

where  X  is  the  real  part  of  one  root  and  the  imaginary  part  (f) 

is  zero 

If  X  is  negative,  Z  approaches  zero  as  time  increases  indefinitely 
and  the  motion  is  said  to  be  overdamped  Like  the  oscillatory  motion,  roots 
which  give  overdamped  motions  will  also  occur  in  pairs.  However,  unlike  the 
complex  conjugate  roots  which  result  in  one  oscillatory  motion,  each  real  coot 
is  a  dinstinct  motion.  Therefore,  it  is  possible  for  an  "n"  degree-of-f re^dom 
system  to  have  "2n"  distinct  motions  if  the  system  is  so  heavily  damped  that 
all  the  roots  to  the  characteristic  equation  are  real. 

There  is  a  third  possible  motion  which  is  a  borderline  case.  If  two 
roots  are  real  and  equal,  the  system  is  said  to  be  critically  damped.  The 
motion  will  be  aperiodic  <ind  both  roots  will  give  the  same  motion. 

The  general  solution  to  the  motion  of  the  system  is  a  linear  combina¬ 
tion  of  all  the  motions  defined  by  the  roots  to  the  characteristic  equation. 
Associated  with  each  root  is  a  mode  shape  which  gives  the  relative  amplitudes 
of  each  degree  of  freedom  when  the  system  is  responding  to  one  particular 
root.  It  is  of  interest  to  establish  these  mode  shapes  so  that  each  stability 
curve  can  be  associated  with  a  definite  motion  of  the  whole  system.  For 
example,  one  mode  shape  may  show  that  the  pitching  motion  of  the  balloon  is 
very  large  compared  to  the  motion  of  the  tether 

3.  STABILITY  ANALYSIS 

Derivations  of  the  equations  of  motion  of  the  tethered  balloon  system 
and  development  of  the  characteristic  equations  for  a  tethered  balloon  system 
approximated  with  a  three-link  tether  are  given  in  Reference  2. 

The  four  linearized  longitudinal  equations  are  Laplace  transformed, 
and  an  eighth  order  characteristic  equation  generated  which  specifies  stability 
characteristics  of  the  system.  In  like  manner,  the  five  linearized  lateral 
equations  can  be  reduced  to  a  tenth  order  equation  which  gives  stability 
information  in  the  lateral  degrees-of-f reedom.  The  roots  of  these  character¬ 
istic  equations  identify  the  natural  frequencies,  damped  frequencies  and 
damping  ratios. 
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Results  of  characteristic  equations  analysis  can  be  displayed  in  the 
form  of  plots  in  a  complex  plane  typically  as  shown  in  Figure  10  of  this 
report.  In  these  plots,  the  abscissa  is  the  real  part  of  the  roots  to  the 
characteristic  equation,  the  ordinate  is  the  imaginary  part.  A  negative  real 
part  means  that  mode  of  oscillation  is  converging  or  stable;  a  positive  real 
part  is  a  diverging  mode.  Only  the  first  and  second  quadrants  are  displayed 
because  the  roots  are  complex  conjugates  which  are  symmetric  to  the  real  axes, 
except  in  the  case  of  overdamped  roots  which  lie  on  the  real  axes.  In  general, 
the  following  information  is  easily  available  for  each  root  directly  from  the 
plots.  The  natural  frequency  is  measured  in  rad/sec  as  the  distance  from  the 
origin  to  the  root;  the  damped  frequency  (rad/sec)  is  the  value  of  the  imaginary 
part  of  the  root;  the  damping  ratio  is  equal  to  the  absolute  value  of  the  real 
part  of  the  root  divided  by  the  natural  frequency. 

4.  DYNAMIC  RESPONSE  ANALYSIS 

The  calculation  of  the  balloon  system  response  to  specific  disturbances 
is  the  subject  of  the  dynamic  response  analysis.  The  most  general  motion  the 
system  can  have  is  a  linear  superposition  of  the  normal  modes. 

Each  aperiodic  or  nonoscillatory  normal  mode  has  one  arbitrary  constant 
(the  initial  value  of  an.  ne  of  the  variables)  associated  with  it;  and  each 
periodic  or  oscillatory  normal  mode  has  two  arbitrary  constants  (the  amplitude 
and  phase  angle  of  any  one  of  the  variables)  associated  with  it.  The  total 
number  of  arbitrary  constants  is  then  equal  to  the  number  of  aperiodic  modes 
plus  twice  the  number  of  periodic  modes;  i.e.  to  the  degree  of  the  character¬ 
istic  equation,  or  the  order  of  the  system.  A  specific  disturbance  will 
excite  the  normal  modes  in  varying  degrees  and  establish  the  values  of  the 
arbitrary  constants. 

The  dynamic  response  of  tethered  balloon  systems  to  various  wind  dis¬ 
turbances  is  obtained  by  integrating  numerically  the  longitudinal  and  lateral 
equations  of  motion  to  produce  a  time  history  of  the  dynamics.  The  start 
conditions,  or  equil  brium  conditions  for  the  dynamic  response  computer  pro¬ 
grams  are  obtained  fi  om  the  linearized  stability  computer  programs  (Reference 
2).  This  approach  to  analysis  has  the  advantage  that  wind  gusts  can  be  pro¬ 
duced  and  the  actual  motion  of  the  system  can  be  observed.  The  major  dis¬ 
advantage  is  that  a  greater  amount  of  computer  time  is  required  when  compared 
to  evaluation  of  stability  by  investigating  the  roots  of  the  characteristic 
equations . 

The  equations  of  motion  for  the  longitudinal  dynamics  of  a  tethered 
balloon  system  were  initially  derived  in  two  forms  (see  Appendix  A):. 

1)  inertia  terms  which  contain  products  of  angular  velocities  are  neglected, 

2)  inertia  terms  which  contain  products  of  angular  velocities  are  included. 

The  concept  of  neglecting  products  of  angular  velocities  is  associated  with 
the  assumption  that  angular  velocities  are  small;'  and  therefore,  products  of 
angular  velocities  are  negligible. 

Numerical  integrations  were  made  with  the  computer  to  determine  the 
effect  of  neglecting  the  inertia  terms  containing  products  of  angular 
velocities..  Although  the  effect  is  obviously  present  (Figures  7  and  8),  the 
overall  difference.1  between  the  results  of  the  two  sets  of  equations  is  small. 
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Figure  7  shows  the  time  history  of  the  four  generalized  coordinates  £l>  £2* 

£3,  and  9.  A  horizontal  gust  of  20  feet  per  second  with  a  duration  of  20 
secords  was  applied.  The  dotted  lines  represent  the  angular  displacement 
when  the  products  of  angular  velocities  are  included  in  the  inertia  terms. 

Figure  8  presents  the  range  and  altitude  of  the  balloon  as  a  function  of  time 
and  also  the  variation  in  tether  tension  at  the  winch  and  at  the  bridle  con¬ 
fluence  point.  Since  this  gust  velocity  is  considered  to  be  high,  it  is 
concluded  that  the  difference  between  the  solid  lines  and  the  dotted  lines 
is  maximum. 

It  was  decided  that  dynamic  simulation  studies  would  be  conducted  with 
a  model  which  neglects  products  of  angular  velocities  for  three  reasons.  First, 
the  equations  containing  products  of  angular  velocities  are  shown  to  give  only 
slightly  different  results.  Second,  it  is  desirable  to  keep  the  dynamic  equa¬ 
tions  compatible  to  the  equations  used  in  the  stability  study  (the  stability 
study  used  linearized  equations).  Third,  it'  is  desirable  to  keep  the  longitu¬ 
dinal  equations  compatible  to  the  l-'teral  equations  (to  derive  the  lateral 
equations  of  motion  containing  products  of  angular  velocities  in  the  inertia 
terms  would  be  a  very  difficult  task  because  of  the  number  of  terms  involved.) 
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SECTION  IV 


DYNAMIC  RESPONSE  OF  TETHERED  BALLOONS 


1.  DESIGN  CONDITIONS  INVESTIGATED 

The  parametric  study  conditions  found  in  Appendix  B,  Tables  B-I  and 
B-II  and  Reference  2  are  the  source  of  information  for  the  results  and 
recommendations  of  this  report.  Each  parametric  dynamic  response  study 
condition  has  a  computer  drawn  plot  of  various  balloon  parameters  versus 
time.  These  plots  are  found  in  Appendix  B. 

The  parameters  chosen  for  investigation  during  this  study  are  placed 
in  three  major  categories: 

(1)  "Balloon  on  Station"  performance  parameters 

(2)  "Structural  Design"  performance  parameters 

(3)  "On  Board  Instrumentation"  Design  parameters 

The  "Balloon  on  Station"  performance  parameters  are  defined  as  balloon 
displacements  from  equilibrium  and,  specifically,  are  changes  in  balloon  range, 
altitude,  pitch  angle,  roll  angle,  yaw  angle,  and  lateral  displacement  induced 
by  a  parameter  variation,  such  as  a  change  in  trim  angle. 

The  "Structural  Design"  performance  parameters  of  interest  are  the 
cable  tension  at  winch  and  bridle,  and  the  balloon  translational  and 
rotational  accelerations. 

The  "On  Board  Instrumentation"  design  parameters  are  defined  as  the 
balloon  angular  displacements  and  rates. 

The  effects  of  gusts,  trim  angle,  altitude,  tail  size,  and  other 
parameters  on  balloon  performance  parameters  are  reviewed  in  detail  later 
in  this  section  of  the  report.  Table  III  lists  the  cases  compared  and  the 
reason  for  comparison  for  the  longitudinal  and  lateral  cases. 

The  operational  winds  used  for  the  investigation,  as  a  function  of 
altitude,  are  listed  in  the  Introduction  to  this  report.  Gusts  used  for 
dynamic  response  are  added  to  the  operational  winds  and  are  applied  to  the 
balloon.  The  gusts  used  in  analysis  are  discussed  next. 

Consider  the  nominal  altitude  condition  of  10,000  feec  where  the 
operational  design  wind  was  67.5  feet  per  second.  Longitudinal  gusts  were 
applied  where  the  horizontal  wind  increased  by  30  percent,  for  a  discrete 
time,  to  a  value  of  1.3  x  67.5  =  87.6  fps.  Gusts  of  2,  10,  and  20  seconds 
and  infinite  time  durations  were  investigated  with  various  rise  times.  A 
case  was  also  investigated  where  the  30  percent  gust  increment  was  applied 
vertically  upward. 
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Table  111 

A.  Comparison  of  Longitudinal  Dynamic  Cases 


Cases 

Balloon 

Compared 

Type 

Reason  for  Comparison 

1,2, 3, 4, 5 

BJ  Nominal 

Gusts  and  Wind  Effects 

3,21,18 

All 

Comparison  of  balloon  Types 

3,6,7 

BJ  Nominal 

Trim  Angle  Effects  * 

18,19, ZO 

VEE  Balloon 

Trim  Angle  Effects  * 

21,22,23 

GAC  Ba1 ’ Don 

Trim  Angle  Effects  * 

11,3,12 

BJ  Nominal 

Tail  Size  Effect 

3,14,15 

BJ  Nominal 

Winching  Effects 

13,3 

BJ  Nominal 

Reduced  Wind  (40%)  Effect 

10,3 

BJ  Nominal 

Amgal  vs  NOLARO 

16,3 

BJ  Nominal 

Effect  of  Payload  on  Underside  of  Balloon 

17,3 

BJ  Nominal 

Winch  at  5000  It,  and  0-0  Ft. 

8,3,9  . 

BJ 

Operational  Altitude  Effects 

Trim  angle 

changed  by  varying  fore  and  aft  location  of  bridle  confluence  point 

B.  Comparison 

of  Lateral  Dynamic  Cases 

Cases 

Balloon 

Compared 

Type 

Reason  for  Comparison 

24,25,26 

BJ  Nominal 

Effects  of  Gust  Ramp  Time 

27,26 

BJ  Nominal 

Effects  of  50%  Increase  in  Gust  Velocity 

40,37,25 

BJ,  VEE, GAC 

Comparison  of  Balloon  Types 

32,25,33 

BJ 

Tail  Size  Effect 

37,38,39 

VEE  Balloon 

Tail  Size  Effect 

35,36 

BJ  Nominal 

Winch  Effects 

34,25 

BJ  Nominal 

Reduced  Wind  Effect  (100%— ►40%) 

31,25 

BJ  Nominal 

Effects  of  Amgal  vs  NOLARO 

29,25,30 

BJ 

Operational  Altitude  Effects 

The  lateral  gusts  applied  to  the  tethered  balloon  systems  to  investi 
gate  dynamic  response  may  be  generally  interpreted  in  two  ways  as  shown  in 
Figure  9  which  depicts  the  wind  in  a  horizontal  plane. 


Condition  A  -  Build-up 
of  wind  from  side 


Condition  B  -  Shift 
in  wind  direction 
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Vw  steady 

- -pi 


VR**Vw 


jvg 
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Longitudinal  axis 
of  balloon 


Vw  steady 


Figure  9  .  Lateral  Wind  Gusts 


In  Condition  A  the  wind  gust  from  the  side  builds  up  with  time.  This 
is  approximately  equivalent  to  the  Condition  B  where  the  wind  remains  constant 
in  magnitude  and  changes  heading  direction  with  time. 

2.  GENERAL  DISCUSSION  OF  DYNAMIC  RESPONSE  OF  TETHERED  BALLOON  SYSTEMS 

The  non-linear  differential  equations  defining  the  motion  of  a  tethered 
balloon  system  and  the  linearized  equations  used  in  studying  the  stability 
characteristics  of  the  balloon  system  are  given  in  Reference  2  together  with 
the  results  of  the  stability  analysis  of  various  balloon  systems.  Appendix  A 
of  this  report  expands  the  equations  of  motion  of  Reference  2  to  include  wind 
gusts.  The  response  of  the  non-linear  balloon  system  when  subjected  to  wind 
gusts  has  been  programmed  for  computer  solution  using  numerical  integration 
methods.  Various  balloon  systems,  subjected  to  different  types  of  gusts, 
have  been  studied  using  this  computer  program  and  th;'  results  are  found  in 
Appendix  B. 

Tha  response  plots  shown  in  Appendix  B  reveal  that  the  balloon  systems 
are  highly  damped.  This  agrees  with  the  stability  analyses  of  the  balloon 
systems  given  in  Reference  2.  In  fact,  the  motion  is  so  highly  damped  that  the 
lower  frequency  response  modes  are,  in  general,  damped  out  within  two  cycles 
or  less  of  motion.  Thus,  the  high  aerodynamic  damping  involved  In  the  balloon 
systems,  makes  it  very  difficult  to  evaluate  the  response  plots  with  regard 
to  phase  relationships  between  coordinates,  damping  characteristics,  and  the 
frequencies  involved  in  the  motion.  It  becomes  evident  after  looking  at  the 
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response  plots  that  only  the  amplitude  and  direction  of  the  response  can  be 
obtained.  A  mathematical  tool  will  be  required  to  obtain  information  regarding 
phase  relationships,  natural  frequencies,  and  damping  characteristics  of  the 
modal  motions  involved.  Such  a  tool  would  be  very  valuable  for  evaluating 
flight  data  from  tethered  balloons. 


What  then,  is  this  mathematical  tool?  The  response  of  a  highly  damped 
multidegree  of  freedom  system  to  a  known  disturbance  can  be  represented 
mathematically  as  the  sum  of  the  system's  modal  responses.  Assuming  the  system 
follows  the  response  of  a  viscously  damped  system,  a  modal  response  definition 
is  given  by- the  following  equation. 


z. 

1 

=  ^c.e 

total  response 

i=N 

ZT 

ii 

N 

represents 

f 

represents 

Wn 

represents 

Wd 

represents 

*i 
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-  ti 
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Z^  of  a  particular  coordinate  is  then  given  by 


(v  ' (i  UH 1 


sin  (  Ud  t  +  0t)  (1) 

i 


the  number  of  frequencies  involved  in  the  response 

the  damping  ratio 

the  natural  frequency 

the  damped  frequency 

the  phase  angle 

an  amplitude  constant 


Each  modal  response  Z^  has  4  unknowns,  since  (j  ^  is  related  to  un 
by  oid  ^  wn  vl  '  s  2  • 

Thus  any  response  curve  which  is  defined  by  a  sufficient  number  of 
points  can  bo  mathematically  represented  by  Equation  (i)  and  solved  for  the 
unknown  constants  using  numerical  curve  fitting  methods.  No  curve  fitting 
program  is  presently  available  to  the  project  for  use  with  Equation  (1). 
however, there  is  a  curve  fitting  program  available  which  uses  a  Fourier 
Transform  and  performs  a  power  spectral  density  (P.S.D.)  analysis  of  the 
Fourier  Transformation.  This  program  is  only  valid  for  analyzing  systems 
with  no  damping.  The  program  still  gives  a  general  overall  view  of  the 
frequencies  involved  in  the  motion  and  it  is  thought  worthwhile  to  analyze 
the  response  of  two  dynamic  cases  and  to  compare  the  frequencies  found  with 
the  frequencies  obtained  from  the  stability  analysis  of  the  same  balloon 
configurations  given  in  Reference  2.  The  dynamic  cases  chosen  for  frequency 
comparison  are  lateral  dynamic  response  Case  28  and  longitudinal  dynamic 
response  Case  4.  The  cases  represent  a  nominal  BJ  balloon  at  10,000  feet. 
Lateral  response  is  for  a  10  second  duration  step  side  gust  and  longitudinal 
response  is  for  a  20  second  duration  step  horizontal  gust. 
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a.  COMPARISON  OF  THE  NON-LINEAR  LATERAL  DYNAMICS  CASE  28,  WITH  THE 
LINEARIZED  STABILITY  CASE  1 

The  small  displacement  motion  of  a  forced  non-linear  system,  in  the 
time  era  after  the  forcing  function  is  removed,  should  result  in  system 
motions  similar  to  the  sum  of  the  mode  shape  times  participation  factor  for 
each  mode  found  in  the  linearized  stability  analysis  of  the  system.  To 
verify  this,  Lateral  Dynamics  Case  28  and  Lateral  Stability  Case  1  (Reference  2) 
were  chosen  for  comparison.  A  complete  description  of  Lateral  Dynamics  Case  28 
is  found  in  Appendix  B,  and  the  Lateral  Stability  Case  1  description  is  found 
in  Reference  2.  For  convenience  the  results  of  Lateral  Stability  Case  1  are 
summarized  in  Table  IV  and  Figure  10  of  this  report.  In  order  to  analyze 
the  coordinate  motions  of  Case  28,  a  power  spectral  density  analysis  with 
Fourier  Transform  was  performed  on  each  of  the  lateral  coordinate  degrees 
of  freedom;  \J/  ,  <£  ,  <7^,  The  input  coordinate  waveforms  and 

PSD  plots  for  each  degree  of  freedom  are  given  in  figures  11  and  12. 

The  PSD  plots  reveal  that  there  are  two  predominate  modes  involved  in 
the  system  motion.  These  are  underdamped  oscillatory  modes  1  and  2  with 
natural  frequencies  and  damping  ratios  of  fn^  =  0.0051  Hz  and  =  0.511 

for  mode  1  and  f  =  0.136  Hz  and  =  0.207  for  mode  2.  Mode  one 

2 

(0.0051  Hz)  is  the  predominate  mode  with  mode  2  (0.136  Hz)  of  smaller  ampli¬ 
tude  superimposed  on  mode  1.  The  lateral  wind  gust  disturbance  in  Case  28 
is  removed  after  T  =  11  seconds,  therefore  all  comparisons  should  be  made 
after  the  first  peak  is  reached  beyond  11  seconds.  It  should  be  noted  that 
mode  2  is  predominately  a  0  or  roll  mode,  and  this  is  evident  from  the  plots 
of  the  coordinates,  since  only  the  0  coordinate  plot  has  evidence  of  a 
0.136  Hz  mode.  On  the  0  plot  mode  2  appears  damped  out  of  the  0  motion  by  the 
time  T  =  37  sec. 


Table  IV.  Summary  of  Lateral  Stability  Case  1 
Normalized  Modes 


Mode 

Frequency 

Damping 

0 

*1 

0 

> 

*3 

Number 

Hz 

Ratio 

Ampl . 

4° 

Ampl , 

4° 

Ampl 

4° 

Ampl. 

in 

Ampl. 

1 

0.0051 

0.511 

1.0 

0.0 

K 

301 

0.342 

83.3 

1WEW 

222 

0.483 

275 

2 

0.1363 

0.237 

1.0 

0.0 

9.12 

252 

0.016 

17.7 

0.051 

215 

0.102 

58.9 

3 

0.000215 

1.0 

1.0 

0.0 

0.0082 

0 

8.86 

180 

2.51 

180 

0.487 

180 

4 

0.0252 

1.0 

1.0 

0.0 

0.927 

0 

1.22 

360 

2.27 

180 

0.957 

360 

5 

0.0578 

1.0 

1.0 

0.0 

1.58 

0 

3.59 

180 

1.54 

360 

0 

0.1295 

1.0 

1.0 

0.0 

23.1 

180 

0.345 

180 

0.147 

360 

0.513 

360 

7 

0.1608 

1.0 

1.0 

0.0 

60.7 

360 

7.34 

180 

6.09 

180 

12.7 

360 

8 

0.4320 

1.0 

1.0 

0.0 

0.988 

360 
.  . . 

0,00007 

180 

0.0006 

360 

0.0051 

180 
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Figure  10.  Lateral  stability  Characteristics  of  Nominal 
BJ  Tethered  Balloon  at  10,000  Feet  Altitude 
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Figure  12.  Lateral  Case  28,  Power  Spectral  Density 


The  aerodynamic  damping  characteristics  of  the  system  appear  viscous 
in  nature  for  the  0.0051  Hz  mode. 

Table  V  lists,  for  Case  28,  the  time  and  amplitude  of  the  various 
coordinates  for  the  first  and  second  peaks  after  the  gust  is  removed. 
Included  in  Table  V  is  the  period,  amplitude  ratio,  and  frequency  as  calcu¬ 
lated  from  the  coordinate  data  of  Case  28. 

The  period  and  amplitude  ratio  for  the  coordinates  i Jf,  0,  and  o ^ 
show  very  good  agreement  with  a  0.0051  Hz  viscously  damped  (  £  =  0.51) 
system.  The  period  between  peaks  of  a  0.0051  Hz  viscously  damped  system 
with  £ =  0.51  is  given  by 


228  sec. 


The  amplitude  ratio  between  successive  peaks  is  given  by 


Xn  +  1 

-2  17  £ 

e 

-2  77  (0.51) 

Xn 

e 

.V 

h.  <2 

.V 

-  0.512 

0.025 


Table  V.  Peak  Amplitudes  for  Lateral  Dynamics  Case  28 


Coordinate 

First  Peak 
After  Gvst 

Second  Peak 

After  Gust 

Period 

T 

Frequency 

1 

Ampl.  Ratio 

Xn  +  1 

Xn 

Time 

Ampl . 

Time 

Ampl  .• 

Sec . 

T 

29 

-3.905 

250 

-0.099 

221 

0.0045 

0.0253 

0 

61 

-0.404 

289 

-0.00973 

228 

0.0044 

0.0241 

26 

-0.475 

197 

0.0943 

171 

0.0059 

0.199 

_ Si _ 

31 

1.33 

224 

0.0630 

193 

0.0052 

0.0474 

76 

-0 . 834 

304 

-0.0207 

228 

0.0044 

0.0248 

b.  COMPARISON  OF  NON  LINEAR  LONGITUDINAL  DYNAMICS,  CASE  4,  WITH  LINEARIZED 
LONGITUDINAL  STABILITY,  CASE  1 

A  description  of  Longitudinal  Dynamics,  Case  4,  is  found  in  Appendix  B 
and  Linearized  Longitudinal  Stability,  Case  1,  is  defined  in  Reference  2.  A 
summary  of  the  Longitudinal  Stability,  Case  1,  is  given  in  Table  VI  and 
Figure  13  of  this  report. 

A  Fourier  and  PSD  analysis  is  performed  on  each  of  the  Longitudinal 
Dynamics,  Case  4,  coordinate  degrees  of  freedom  0,  |^,  ?2>  ant*  £3- 

The  input  coordinate  waveform  and  PSD  plot  for  each  degree  of  freedom  is 
given  in  Figures  14  and  15  respectively. 

The  need  of  a  mathematical  tool  for  analyzing  a  highly  damped  response 
is  very  evident  in  Longitudinal  Dynamics,  Case  4.  The  response  of  every 
coordinate  is  damped  out,  or  nearly  so,  in  less  than  one  cycle.  The  Fourier 
and  PSD  non-damped  harmonic  ana lyis  of  this  motion  is  not  too  meaningful. 
However,  it  does  indicate  that  the  motion  is  comprised  of  frequencies  below 
0.05  Hz  with  the  most  power  below  0.01  Hz. 


Table  VI.  Summary  of  Longitudinal  Stability  Case  1 

(Reference  2) 


Norma 1 i z ed  Modes 


Mode 

Number 

Frequency 

Hz 

Damping 

Ratio 

e 

« 

_ _ 

? 

I-] 

Ampl.' 

Ampl, 

EM 

Ampl. 

EM 

Ampl. 

WOW 

1 

0.00554 

-0.037 

1.0 

0.0 

2.14 

281. 

1.27 

238. 

1.28 

194.0 

mm 

0.0728 

0.653 

1.0 

mwm 

0.131 

197. 

0.080 

353. 

0.058 

39.7 

0.0986 

0.51 

1.0 

0.0 

0.0747 

9.6 

0.214 

188. 

0.145 

0.0156 

1.0 

1.0 

0.0 

0.888 

0.0 

0.800 

mm 

0.363 

ieo. 

0.1715 

1.0 

1.0 

0.0 

0.00033 

■a 

0.0123 

0.0 

0.0066 

180. 

t.  YAW-ROLL  COUPLING  IN  LATERAL  DYNAMIC  CASES 

Yaw-roll  coupling  for  BJ  and  GAC  single  hull  balloons  subjected  to 
lateral  wind  is  evident  when  the  \J/  response  plots  are  reviewed.  The  Vee 
balloon  appears  to  have  negligible  coupling.  As  an  example,  observe  the 
motion  of  coordinate  xfy  (Reference  Appendix  B)  for  Lateral  Case  28.  One 
would  think  that,  when  a  balloon  is  in  equilibrium  with  a  steady  head  wind 
and  then  is  subjected  to  a  side  gust  from  the  left,  the  balloon  would  yaw 
into  the  resultant  wind.  This  would  be  a  negative  yaw  angle,  however, 
Lateral  Case  28  shows  the  balloons  first  yaw  motion  is  positive  or  away  from 
the  resultant  wind  for  a  time  of  approximately  7  seconds  and  to  =  +4° 
before  the  balloon  starts  yawing  toward  the  resultant  wind. 
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Figure  13,  loxgtudlnal  Stability  Characteristics  of  Ifarioal 
BJ  Tethered  Balloon  at  10,090  Feet  Altitude 


w 


To  explain  this  motion  refer  to  Equation  86  in  Appendix  A.  Neg¬ 
lecting  link  acceleration  terms  results  in  the  following  equation. 

V  -  |MS  [-*J.  e  33  +  *10.  f  23]  [-  y]  -  *YB  59  +  *YZB  09  «•}  9 

+  |Ks[-V£  33  "  »kn«2|  +[«l  -V‘13] 

+  K  Sjm  *  uj  +  C2  9  j"lXB  S2  0  +  *ZB  C  d] 

+  S2  0  IYB  -  S29  Cd  Ijz,|» 

#  *  •  • 

Substituting  a  and  b  for  the  coefficients  of  0  and  \p  results  in  the 
following  equation. 

•  • 

F^  =  a  0  +  bi|j 

The  equation  is  solved  for 


b 

where  F  ^  is 

the  total  externally  applied  balloon  y 

omen  t . 

Case  28 

parameters  have  the  following  values  at 

T  =  5.5  seconds 

(from  digital 

printout  of  Case  28). 

II 

pH 

-40039  ft- lb 

0  =  3.060  °/sec2 

2 

ip  =  1.207°/sec 

*XB  = 

116,000  Slug  ft2 

0  =  2.373  7sec 

Ip  -  1 . 989°/ sec 

IYB  = 

26,000  Slug  Ft2 

0  =  0.531° 

lp  =  0.426° 

XZB 

119,000  Slug  Ft2 

IYZB 

12,500  Slug  Ft2 

28 


i 


t 


i 

: 

7 
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Inspection  of  the  Jjj  coefficient,  b,  reveals  it  to  be  positive.  Since 
F<J>  *  -40039  ft- lb  and  will  tend  to  turn  the  balloon  into  the  resultant  wind, 
the  only  term  left  to  cause  i|>  to  be  positive  is  the  roll  coupling  term  a<£ 

In  conclusion  the  initial  ij>  motion  away  from  the  resultant  wind  is  caused 
by  roll  (0)  inertia  coupling. 

3.  EFFECTS  OF  DESIGN  PARAMETERS  ON  TETHERED  BALLOON  DESIGN 
a.  Balloon  on  Station  Performance  Parameters 

(1)  General.  As  listed  in  Part  1  of  this  section  a  number  of  tethered 
balloon  design  parameters  and  flight  conditions  were  varied  to  establish  their 
influence  on  a  tethered  balloon  design.  This  subsection  summarizes  some  of 
the  major  effects  on  tethered  balloon  dynamic  behavior  as  established  from 
the  dynamic  response  calculations  presented  in  Appendix  B  and  an  evaluation 
of  these  data. 


(2)  Wind  Gust  Effects,  "he  response  data  in  Table  VII  shows  the  effect 
of  increasing  the  duration  of  a  horizontal  gust  on  the  longitudinal  motions 
of  the  nominal  BJ  balloon  design  for  10,000  foot  altitude.  It  is  apparent 
from  an  examination  of  these  data  that  the  longer  duration  of  the  gust  the 
greater  the  angle  of  attack  of  the  balloon  becomes.  A  new  equilibrium  con¬ 
dition  is  reached  when  the  'ust  is  left  on  indefinitely.  As  a  consequence 
of  the  increased  angle  of  attack  and  wind  velocity  the  aerodynamic  lift 
increases  and  a  somewhat  more  favorable  aerodynamic  lift  to  drag  ratio  is 
obtained.  This  serves  to  increase  the  altitude  of  the  balloon  and  reduce 
the  downwind  displacement  of  the  balloon.  A  substantial  increase  in  tether 
tension  is  also  evident. 


s 

f 
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A  nominal  BJ  balloon  operating  at  10,000  feet  is  subjected  to  a 
20.2  ft/sec,  20  second  deration,  horizontal  gust  in  Longitudinal  Case  4  and 
to  a  20.2  ft/sec,  20  second  duration,  vertical  gust  in  Longitudinal  Case  4A. 
Table  VII  also  summarizes  the  results  of  these  two  cases. 


Table  VII.  Wind  Gust  Duration  Effects  on  BJ  Nominal  Balloon 

in  Longitudinal  Plane 


Altitude  “  10,000  Ft.  Gust  Velocity  ■  20.2  Ft/Sec 

Payload  Weight  •  1,000  Lb.  Gust  Duration  *  tr 

Trim  Angle  -8.5°  0 


Balloon  Displacements 

From  Equilibrium 

Cable  Tensions  At 
Equilibrium 

Maximum  Cable 

Tensions 

t 

A 

A 

A 

Case 

G 

Range 

Altitude 

Pitch 

Winch 

Bridle 

Winch 

Brid le 

No., 

Down/ Up 

Sec . 

Ft.; 

Ft.- 

Deg.- 

Lb. 

Lb. 

Lb, 

Lb. 

1 

OO 

-65/800 

690 

5.7 

1650 

2210 

3630 

A200 

2 

2 

-12/16 

15 

1.9 

1650 

2210 

2250 

2830 

3 

10 

-52/73 

82 

2.9 

1650 

2210 

2650 

3210 

4 

20 

-t>8/]46 

121 

A. 7 

1650 

2210 

3180 

3740 

4A 

20 

0/A51 

3A5 

_ LiiJ 

1650 

2210 

3582 

394  7 

1 

\ 
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Specific  comments  for  this  comparison  are: 

(1)  The  vertical  gust  displaced  the  balloon  up  range  approxi¬ 
mately  three  times  the  horizontal  gust  up  range  displacement. 

(2)  The  vertical  gust  caused  the  balloon  to  rise  three  times  higher 
than  the  horizontal  gust. 

(3)  The  pitch  angle  change  from  a  vertical  gust  is  only  1/4  the 

pitch  angle  change  due  to  a  horizontal  gust  of  equal  magnitude. 

(4)  Maximum  cable  tension  increased  by  400  pounds  due  to  the 
vertical  gust.  It  is  noted  that  equilibrium  cable  tension  is 
slightly  different  from  that  listed  in  Table  I.  This  is  a 
result  of  different  balloon  trim  angles  and  number  of  cable 
links  used  in  the  two  analyses. 

The  effects  of  gust  ramp  rise  time  on  the  response  of  a  nominal  BJ 
balloon  is  given  in  Table  VIII.  These  cases  represent  BJ  nominal  balloons 
at  10,000  feet  chat  are  subjected  to  a  steady  lateral  gust  with  ramp  rise 
times  of  1,  3  and  5  seconds. 

Table  VIII  reveals  that  gust  linear  ramp  rise  times  of  1,  3,  and  5 

seconds  have  little  effect  on  balloon  response,  except  for  balloon  roll 

where  the  1  second  rise  time  is  1.8  times  larger  than  the  5  second  rise 
time  response. 


Table  VIII.  Wind  Gust  Ramp  Rise  Time  Effects  on  BJ  Nominal  Balloon 

in  Lateral  Plane 


Cable  Tension 

At  EquI librium 

Ramp 

Gust 

from  Initial  Equillt 

rlum 

Ten 

lion 

Case 

Rite 

Velocity 

Yaw 

Roll 

Ut. 

At 

At 

At 

At 

No. 

Time 

Ft/Sec 

+/- 

+/- 

Displ 

Winch 

Bridle 

Winch 

Bridle 

Deg. 

Deg. 

Ft.' 

Lb. 

Lb. 

Lb. 

Lb. 

24 

1  sec 

12 

4  /-13.6 

3.9/-0.4 

460 

1650 

2210 

1686 

2258 

25 

3  -ec 

12 

4.3/-13. 8 

3.1/-0.35 

460 

1650 

2210 

1680 

2250 

26 

5  se^ 

12 

4  /- 13.6 

2.2/-0.4 

460 

1650 

2210 

1675 

2244 

The  effects  of  a  50  percent  increase  in  lateral  gust  velocity  on  the 
response  of  a  nominal  BJ  balloon  is  given  in  Table  IX.  The  balloon  is  opera¬ 
ting  at  10,000  feet  in  a  steady  head  wind  when  subjected  to  steady  lateral 
gusts  of  12  ft/sec  and  18  ft/sec.  Each  gust  has  a  ramp  rise  time  of  five 
seconds.  Table  IX  reveals  that  a  50  percent  increase  in  steady  lateral  gust 
velocity  results  in  approximately  a  50  percent  increase  in  balloon  lateral 
displacement,  roll,  and  yaw  response. 


30 


_-y Z?  - 


Table  IX.  Effect  of  Increase  in  Gust  Velocity  on  BJ  Nominal 
Balloon  in  Lateral  Plar' 


Balloon  Displacements 

Cable  Tension 

Maximum  Cible 

from  Initial  Equilibrium 

At  Equilibrium 

Tension 

Yaw 

Roll 

Lat. 

At 

At 

At 

At 

Case  No. 

+/- 

+/- 

Displ. 

Winch 

Bridle 

Winch 

Bridle 

Deg. 

Deg. 

Lb. 

Lb. 

Lb. 

Lb. 

26 

4.0/-13.6 

2.2/-0.4 

+60 

1650 

2210 

1675 

2244 

27 

6.0/-20.0 

3.2/-0.6 

710 

1650 

2210 

1708 

2285 

(3)  Comparison  of  Balloon  Types.  The  effects  of  trim  angle  on  the 
longitudinal  response  of  BJ  nominal,  VEE,  and  GAC  single  hull  balloons  is 
given  in  Table  X.  The  balloons,  in  equilibrium  with  a  steady  head  wind  at 
10,000  feet,  are  subjected  to  a  horizontal  longitudinal  gust  of  20.2  feet/ 
second  for  ten  seconds.  Table  X  reveals  that: 

(1)  the  down  range  displacements  from  equilibrium  decreases 
with  an  increase  in  trim  angle  of  attack 

(2)  the  GAC  single  hull  balloon  has  the  smallest  displacement 
from  equilibrium  of  the  three  balloons 

(3)  the  BJ  nomi  ,.1  balloon  has  the  smallest  cable  tension  of  the 
three  balloon  types 

(4)  cable  tension  increases  with  an  increase  in  trim  angle 

(5)  changes  in  trim  angle  have  a  negligible  effect  on  pitch 
response  to  gusts. 

The  effects  of  lateral  steady  wind  gusts  on  the  lateral  response  of 
BJ  nominal,  VEE,  and  GAC  single  hull  balloons  are  given  in  Table  XI.  The 
balloons  are  in  equilibrium  at  10,000  feet  with  a  steady  head  wind  and  then 
subjected  to  a  steady  lateral  horizontal  gust  of  12  ft/sec  with  a  three 
second  ramp  rise  time.  Table  XI  reveals  that: 

(1)  the  response  of  the  three  balloons  is  approximately  the  same 
in  both  yaw  and  roll  response 

(2)  the  VEE  balloon  has  1/3  the  lateral  displacement  of  the  other 
balloons 

(3)  the  VEE  balloon  has  2.5  times  more  tether  load  than  the  BJ 
balloon  as  a  result  of  the  larger  aerodynamic  lift  of  this 
balloon. 


Table  X.  Effect  of  Trim  Angle  of  Attack  on  Longitudinal  Response  of 

Three  Balloon  Types 


Balloon  Displacements 
from  Equilibrium 

Cable  Tenalone 

At  £qui librium 

laxlanoi  Cable 
Tensions 

Type 

Case 

Trim 

D  A 

Range 

& 

Pitch 

At 

At 

At 

At 

Balloon 

No. 

Angle 

Winch 

Bridle 

Winch 

Bridle 

Ft.; 

Ft.- 

Degree 

Lb., 

lb. 

lb. 

lb. 

BJ 

6 

Q 

H 

65 

2.8 

1350 

1880 

2070 

2610 

Norn. 

3 

■ 9 

Bl 

82 

2.9 

1650 

2210 

2650 

3210 

7 

10.5 

-50/75 

55 

2.9 

1950 

2525 

3220 

3800 

n 

mg 

-35/135 

80 

4.2 

BIS 

9500 

7500 

VEE 

El 

la 

-22/119 

60 

3.7 

m'M 

5802 

10100 

20 

9.1 

-19/105 

A0 

3.3 

6000 

7200 

12900 

14000 

CAC 

22 

6.0 

m 

MM 

mmm 

EB 

3070 

3150 

4000 

Single 

21 

7.9 

mM 

■SB 

I 

3596 

3940 

4830 

Hu  11 

23 

10.0 

-90/92 

wm 

mm 

mm 

4150 

4B50 

5770 

Table  XI.  Effect  of  Lateral  Wind  Gusts  on  the  Lateral  Response 
of  the  Three  Balloon  Types 


Balloon  Displacements 
from  Equilibrium 

Cable  Tensions 

At  Equilibrium 

Maximum  Cable 
Tensions 

Type 

Case 

Yaw 

Roll 

Lat . 

At 

At 

At 

At 

Balloon 

No. 

+/- 

+/- 

DiSple' 

Winch 

Bridle 

Winch 

Bridle 

Deg. 

Deg.. 

Ft, 

lb. 

lb. 

Lb. 

Lb, 

BJ 

25 

4.3/-13.8 

3.1/-0.35 

460 

1650 

2210 

1680 

2250 

VEE 

37 

0.5/- 15 . 5 

3. 1/-0.4 

178 

4620 

5802 

4683 

5B22 

GAC 

40 

5.5/- 14 . 5 

2. 7/-0.5 

510 

2696 

3596 

2734 

3645 

r,  y*m*~*' 


y 

i 

lr 

K 


(4)  Effect  of  Tail  Size.  The  effect  of  tail  size  on  the  response  of 
nominal  BJ  balloons  subjected  to  longitudinal  gusts  is  given  in  Table  XIX. 

The  effects  of  tail  size  on  the  response  of  BJ  nominal  and  VEE  balloons 
subjected  to  lateral  gusts  are  given  in  Table  XIII.  The  longitudinal  cases 
of  Table  XII  are  in  equilibrium  with  a  steady  headwind  and  then  subjected 

to  20.2  ft/sec  step  gust  of  10-second  duration.  The  tail  sizes  investigated 
were  81%,  100%  and  144%  of  nominal  tail  area.  Table  XII  reveals  that: 

(1)  both  up  and  down  range  displacements  decrease  with 
increasing  tail  size 

(2)  pitch  angle  response  decreases  with  increasing  tail  size 

(3)  the  maximum  tether  load  at  bridle  and  winch  decrease  with 
increasing  tail  size. 

The  lateral  cases  of  Table  XIII  are  in  equilibrium  with  a  steady  head 
wind  and  then  subjected  to  a  lateral  steady  gust  of  12  ft/sec  with  a  ramp 
rise  time  of  three  seconds.  Table  XIII  reveals  that: 

(1)  the  BJ  balloon  has  larger  lateral  displacement  response 
and  less  roll  response  with  increasing  tail  size 

(2)  the  VEE  balloon  has  smaller  response  in  yaw,  roll,  and 
lateral  displacement  with  increasing  tail  size 

(3)  tether  loads  for  both  types  of  balloons  increase  with  an 
increase  in  tail  size. 

(5)  Effects  on  Balloon  Behavior  of  Intermediate  Altitudes.  The  effects 
of  intermediate  altitude  operation  on  a  nominal  BJ  balloon  designed  for 
10,000  feet  are  given  in  Table  XIV.  The  altitudes  investigated  are  1,000, 

4,000  and  10,000  feet.  The  balloon  is  in  equilibrium  with  a  steady  head 
wind  and  then  is  subjected  to  longitudinal  gusts  of  10-second  duration  with 
varying  velocities  of  20.2  ft/sec,  13.6  ft/sec,  and  7.1  ft/sec  for  longitudinal 
dynamics  cases  3,  14,  and  15,  respectively.  The  wind  gust  velocities  are 
reduced  proportional  to  the  steady  wind  velocity  reductions  with  altitude. 

Table  XIV  reveals  the  longitudinal  responses  all  decrease  with  a  decrease  in 
altitude  andthe  maximum  tether  load  at  bridge  is  proportional  to  the  gust 
velocity. 

The  effects  of  intermediate  altitude  operation  on  the  lateral  response 
of  a  BJ  balloon  designed  for  10,000  feet  are  given  in  Table  XV.  The  balloon 
is  in  equilibrium  with  a  steady  head  wind  and  then  is  subjected  to  lateral 
steady  gusts  of  12  ft/sec,  8  ft/sec,  and  1  ft/sec  corresponding  to  altitudes 
of  10,000,  4,000,  and  1,000  feet.  Table  XV  reveals  that:. 

(1)  the  lateral  responses  of  the  nominal  BJ  balloon  subjected  to 
lateral  steady  gusts  increase  with  decreasing  altitude 

(2)  a  review  of  the  computer  plots  of  Case  36  (1000  feet)  in 
Appendix  B  indicates  the  response  is  unstable. 

A  detailed  discussion  of  Lateral  Case  36  follows. 
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Table  XII.  Longitudinal  Response  of  BJ  Balloon  With  T’il 

Size  Variation 


Ballbnn  Displacements 

From  Equilibrium 

Cable  Tensions 

At  Equilibrium 

Maximum  Cable 
Tensions 

Trim 

Tall 

Case 

Range 

Alt.. 

Pitch 

At 

At 

At 

At 

Angle 

Size 

No.. 

Down/Up 

Winch 

Bridle 

Winch 

Bridle 

Deg. 

Ft. 

Ft. 

Deg. 

Lb. 

Lb. 

Lb. 

Lb. 

11.0 

1 

11 

-61/88 

70 

4.6 

1730 

1880 

3000 

3550 

8.5 

3 

-52/73 

82 

2.9 

1650 

2210 

2650 

3210 

7.4 

144T 

12 

-50/55 

40 

1.3 

1750 

2320 

2600 

3170 

Table  XIII.  Lateral  Response  of  BJ  and  7EE  Balloons  With 
Tail  Size  Variation 


Balloon  Displacements 
from  Equilibrium 

Cable  Tensions 
at  Equilibrium 

Maximum  Cable 
Tensions 

Type 

Tail 

Trim 

Case 

Yaw 

Roll 

Lat. 

At 

At 

Balloon 

Size 

Angle 

No., 

+/- 

+/- 

Displ ., 

Winch 

Bridle 

Deg. 

Deg.. 

Ft, 

Lb. 

Lb. 

8  IT 

10.98" 

32 

3.4/- 13. 5 

11/- 1 

420 

1700 

mm 

B 1 

2400 

BJ 

1007. 

8.5" 

25 

4.3/- 13. 8 

3, 1/-0.35 

460 

1650 

ISII&hI 

m 

22bO 

144% 

_ 

7.37* 

33 

6.1/-13.3 

2.3/-0.6 

480 

2383 

2901 

2435 

2972 

m 

6.96* 

37 

0.5/-15 

3.1/-0.4 

177 

4650 

4605 

4683 

4855 

VEE 

7.23° 

38 

0  /-U 

1.9/-0.2 

78 

4763 

| 

4771 

5909 

II 

7.5° 

39 

0  /-ll 

0.7/-0.1 

72 

4874 

■  I 

4883 

6021 

Table  XIV.  Effect  of  Intermediate  Altitudes  on  Longitudinal 

Response  of  the  BJ  Balloon 


Balloon  Displacements 
from  Equilibrium 

— 

Cable  Tensions 
at  Equilibrium 

Maximum  Cable  1 
Tension  | 

Cese 

Altitude 

Gust 

Renge 

Alt. 

Pitch 

At 

At 

At 

At 

No. 

Velocity 

Down /Up 

Winch 

Bridle 

Winch 

Bridle 

Feet 

Ft/Sec 

Ft.. 

Ft. 

Deg. 

Lb.. 

Lb. 

lb.. 

Lb. 

3 

10,000 

20.2 

-52/73 

82 

2.9 

165(. 

2210 

2650 

3210 

14 

4,000 

13.6 

-30/0 

-10 

0.6 

1420 

1650 

2040 

2260 

15 

1,000 

7.1 

-25/10 

-  3 

0.45 

980 

1020 

1110 

1160 

Table  XV. 


Effect 


of  Intermediate  Altitudes  on  Lateral  Response 
of  the  BJ  Balloon 


Balloon  Disp^cements 
from  Enui librium 

Cable  Tensions 
at  Equilibrium 

Maximum  Cable 
Tension* 

Case 

Gust 

Alt. 

Yaw 

Roll 

Let. 

At 

At 

At 

At 

No. 

Velocity 

>/- 

+/- 

Displ. 

Winch 

Bridle 

Winch 

Bridle 

Ft/Sec 

Ft.. 

Deg.. 

Deg. 

Ft. 

_  Lb'  ,_J 

Lb.. 

Lb. 

25 

12.0 

10,000 

4.3/-13.8 

3.1/-0.35 

460 

1650 

2210 

1680 

2250 

35 

8.0 

4,000 

13. 1/-27.0 

2.1/-1.26 

490 

1422 

1645 

1449 

1667 

36 

1.0 

1,000 

52./-63.3 

10.2/-5.5 

570 

1031 

1961 

1988 

Case  36  represents  a  nominal  BJ  Balloon  designed  for  a  10,000-foot 
operational  altitude  and  operating  at  1,000  feet  altitude.  The  balloon  is  in 
equilibrium  with  a  steady  headwind  and  then  is  subjected  to  a  1  foot/second 
steady  lateral  wind. 

A  lateral  stability  analysis  of  this  condition  is  found  in  Reference  2 
Run  No.  22.  Reference  2  (Figure  42)  reveals  an  unstable  mode  with  a  frequency 
of  0.011  Hz.  The  shape  of  this  mode  is  also  found  in  Reference  2  (Figure  11) 
and  has  the  following  relationships  between  amplitudes  and  phase  angles  when 
the  mode  is  normalized  on  the  balloon  yaw  angle. 


Coordinate 

Relative  ^ 
Amplitude 

Phase  Angle 
Degree 

Balloon  Yaw 

1.000 

0.0 

Balloon  Roll 

0 

0.253 

101.3 

Link  1  Yaw 

0.340 

82.4 

Link  2  Yaw 

°  l 

0.332 

81.9 

Link  3  Yaw 

®  3 

0.324 

80.9 

Data  is  taken  from  computer  digital  printout  of 
Run  22  (this  data  not  included  in  Reference  2). 

The  stability  analyses  also  reveals  this  mode  has  a  rise  time  to 
double  amplitude  of  86.2  seconds. 

Turning  to  the  computer  plots  for  lateral  case  36  (see  Appendix  B) , 
the  following  observations  are  made: 

(1)  Tether  links  are  in  phase  and  of  approximately  equal 
amplitude . 

(2)  The  balloon  roll  angle  is  approximately  90°  out  of  phase 
with  respect  to  the  balloon  yaw  angle. 

(3)  The  three  links  are  approximately  83  degrees  out  of  phase 
with  respect  to  the  balloon  yaw  angle. 

(4)  The  rise  time  to  double  amplitude  for  the  balloon  yaw  angle 
is  approximately  85  seconds. 

(5)  The  oscillation  has  a  period  of  93  seconds  (0.0109  Hr,). 

(6)  Observations  (1)  through  (5)  indicate  that  lateral  case  36 
is  predominantly  the  unstable  0.011  Hz  mode  found  in  the 
stability  analysis  run  22. 

(7)  It  should  be  noted  also  that  the  oscillations  appear  to  be 
diverging,  as  would  be  expected,  since  the  0.011  Hz  mode  is 
predominate . 
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In  conclusion  lateral  case  36  indicates  a  nominal  BJ  Balloon  designed 
for  a  10,000-foot  operating  altitude  will  have  lateral  instability  problems, 
if  it  is  operated  at  1,000  feet.  Lateral  cases  35  and  37  fly  the  same 
balloon  at  4,000  and  10,000  feet  and  reveal  no  diverging  motion.  GAC  has 
investigated  this  type  of  instability  problems  in  Reference  4.  In  this 
Reference  the  formulation  of  a  method  of  analysis  to  predict  the  instability 
of  short  tethered  balloons  exposed  to  wind  was  derived  and  experimental  tests 
were  performed  on  towed  models  to  verify  the  analysis.  As  a  result  of  this 
work,  it  can  be  stated  that  lateral  instability  occurs  when  the  inverted 
pendulum  frequency  of  the  balloon  on  the  tether  coincides  with  the  balloon 
yaw  frequency.  Since  the  pendulum  frequency  is  inversely  proportional  to 
the  square  root  of  the  tether  length,  it  is  obvious  that  for  a  certain  wind 
velocity  a  critical  tether  length  exists.  It  should  also  be  mentioned  that 
Case  36  is  not  a  winching  condition,  since  the  cable  v?as  a  fixed  length  in 
the  analysis.  However,  it  does  suggest  that  a  possible  stability  problem 
will  occur  when  the  balloon  is  winched  through  the  1,000  foot  altitude  level. 

(6)  Effects  of  Reduced  Winds  (40  Percent  of  Nominal  Winds) .  The  effects 
of  reducing  the  steady  equilibrium  wind  and  gusts  to  40  percent  of  nominal 
on  response  in  the  longitudinal  and  lateral  planes  are  summarized  in  Table  XVI. 
In  comparing  the  effects  of  a  40  percent  reduction  in  winds.it  should  be 
noted  that  the  equilibrium  trim  angle  decreased  from  8.5  degrees  at  100  per¬ 
cent  of  nominal  wind  to  2.96  degrees  at  40  percent  of  nominal  wind.  The  yaw 
response  reveals  no  change  in  response,  while  the  other  responses  are  reduced 
as  expected. 


Table  XVI,  Effect  of  Reducing  Winds  on  Longitudinal  and  Lateral  Response 

a)  Longitudinal  Response 


Balloon  Displacements 
from  Eouilibrium 

Cable  Tensions 
at  Equilibrium 

Maximum  Cable 
Tensions 

Case 

Wind 

Trim 

Range 

Ait.. 

Pitch 

At 

A 

At 

At 

No.' 

And 

Angle 

Down/Up 

Winch 

Bridle 

Winch 

Bridle 

Ousts 

Deg. 

Ft. 

Ft. 

Deg. 

Lb.- 

Lb.. 

Lb. 

3 

Nominal 

6.3° 

-52/73 

62 

2.9 

1650 

2210 

2650 

3210 

13 

40%  Nom. 

2.96J 

-38/0.0 

0.63 

315 

- 

875 
...  ■ 

355 

910 

b)  Lateral  Response 


- . - . — • — ■ - - - 

Balloon  Displacements 

Cable  Tension 

Maximum  Cable 

from  Initial  Equilibrium 

at  Equilibrium 

Tension 

Case 

Wind 

Trim 

Yaw 

Roll 

I.at. 

At 

At 

At 

At 

No. 

And 

Angle 

+/- 

+/-- 

Displ . 

Winch 

Bridle 

Winch 

Bridle 

Gusts 

Deg, 

Deg. 

Deg. 

Ft. 

Lb. 

I.b.. 

Lb. 

Lb. 

25 

Nomina l 

8.5 

4 . 3  /  -  1  5  a 

3. 1/-0.35 

1 

460 

1650 

2210 

1680 

2250 

34 

40^.  Nom . 

_ 

2.96 

4  9/-13.8 

0 . 91 / -0.07 

322 

31b  5 

_ 

871.6 

321.2 

874.7 

37 


(7)  Effects  of  Nolaro  versus  Amgal  Tethers.  The  effects  of  Nolaro 
and  Amgal  tethers  on  the  longitudinal  response  of  BJ  type  balloons  are  given 
in  Table  XVII.  The  BJ  type  balloons  are  at  10,000  feet  in  equilibrium  with 
a  steady  head  wind  and  then  are  subjected  to  a  20.2  ft/sec  horizontal 
longitudinal  gust  of  10  seconds  duration.  The  balloon  using  the  Nolaro  tether 
has  a  volume  of  60,000  cu.  ft.,  while  the  balloon  using  the  Amgal  tether  has 
a  volume  of  75,000  cu.  ft.  The  larger  balloon  is  required  because  the  Amgal 
tether  is  heavier  than  the  Nolaro  tether.  The  trim  angle  for  both  cases  is 
8.5  .  Table  XVII  reveals  the  Nolaro  system  has  slightly  more  response  and 
smaller  tether  loads  than  the  Amgal  system. 


Table  XVIL  Effect  of  Nolaro  and  Angal  Tethers  on  Longitudinal  Response 


— 

Balloon  Displacements 
From  Equilibrium 

— 

Cable  Tensions 

At  Equilibrium 

ivia^ 

Tensions 

Cable 

Case 

Range 

Up/Down 

Alt 

Pitch 

Winch 

Bridle 

Winch 

Bridle 

Type 

No. 

Ft 

Ft. 

Deg. 

Lb 

i_b 

Lb 

Lb 

Nolaro 

WM 

-52/73 

82 

2.9 

1650 

2210 

mm 

mm 

Amgal 

■9 

-60/50 

28 

2.  1 

1880 

2870 

EH 

KH 

The  effects  of  Nolaro  and  Amgal  tethers  on  the  lateral  response  of 
60  K.  Ft^  and  75K.  F^  BJ  type  balloon  sv .terns  are  given  in  Table  XVIII.  The 
baHoon  systems  are  in  equilibrium  with  a  steady  head  wind  and  then  subjected 
to  a  steady  lateral  horizontal  gust  of  12  ft/sec  with  a  ramp  rise  time  of  three 
seconds.  Table  XVIII  reveals  that 

(1)  Yaw  and  roll  response  of  the  Nolaro  system  are  slightly 
larger  than  the  response  of  the  Amgal  system 

(2)  The  lateral  displacement  of  the  Amgal  system  is  1.2  times 
the  response  of  the  Nolaro  system 

(3)  The  Amgal  system  has  larger  tether  loads  as  expected  since 
the  Amgal  tether  is  heavier. 


Table  XVIII.  Effect  of  Nolaro  and  Amgal  Tethers  on  Lateral  Response 


Balloon  Displacements 
From  Initial 

Equilibrium 

Cable  Tension 
At  Equilibrium 

Max  Cable 
Tension 

Trim 

Yaw 

Roll 

Lat 

At 

At 

At 

At 

Case 

Cable 

Angle 

+/- 

+/- 

Displ. 

Winch 

Bridle 

Winch 

Bridle 

No. 

Type 

Deg. 

Deg. 

Deg. 

Ft. 

Lb 

Lb 

Lb 

Lb 

25 

Nolaro 

8.  5 

4.  3/-14 

3.  l/-.  35 

E39 

fM 

2210 

1680 

2250 

31 

Amgal 

8.  3 

4  /- 13 

2.  1/  -..4 

H 

B 

2872 

1905 

2905 

38 
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(8)  Effects  of  Payload  Location.  The  effects  of  payload  location  on 
the  longitudinal  response  of  a  nominal  BJ  balloon  at  10,000  ft  are  given  in 
Table  XIX.  Case  three  has  the  1,000  pound  payload  located  on  the  confluence 
point  of  the  bridle  while  case  sixteen  has  the  payload  mounted  on  the  underside 
of  the  balloon.  The  balloons  are  in  equilibrium  with  a  steady  head  wind  and 
then  subjected  to  a  head  gust  of  20.2  ft/sec  with  a  duration  of  10  seconds. 
Table  XIX  reveals  the  configuration  with  the  payload  at  the  bridle  confluence 
point  (case  3)  has  approximately  40  percent  more  longitudinal  response  than 
the  payload  mounted  on  the  balloon  bottom  configuration.  The  tether  leads 
are  approximately  the  same  for  both  configurations. 


Table  XIX.  Effect  of  Payload  Location  on  Longitudinal  Response 


Balloon  Displacements 
From  Equilibrium 

L. _ 

Cable  Tension 
At  Equilibrium 

Max  Cable 
Tension 

Case 

Payload 

Range 

Down/Up 

At 

W  inch 

At 

Bridle 

At 

Winch 

At 

Bridle 

No. 

Location 

Ft. 

Lb. 

Lb. 

Lb. 

_ 

Lb. 

_ i 

3 

Confl. 

-52  /73 

82 

2.9 

mm 

I 

mm 

W 

16 

Underside 

-49  /56 

48 

2.  1 

B 

B 

(4)  Effect  of  Winch  Altitude  Location.  The  effect  of  winch  altitude 
location  on  the  longitudinal  response  of  a  nominal  BJ  balloon  remaining  at 
10,000  feet  MSL  is  given  in  Table  XX.  The  winch  is  located  at  0.0  ft  and 

5,000  ft  for  cases  3  and  17  respectively.  The  balloons  are  in  equilibrium 

with  a  steady  head  wind  and  subjected  to  a  20.2  ft/sec  gust  cf  10  seconds 

duration.  Table  XX  reveals  that, 

(1)  the  range  and  altitude  responses  are  reduced  when  the  winch 
is  at  5 ,000  feet  , 

(2)  the  pitching  response  is  slightly  larger  with  the  winch  at 
5,000  feet, 

(3)  the  tether  load  is  maximum  when  the  winch  is  ct  5,000  feet. 


Table  XX.  Effect  of  Winch  Altitude  Location  on  Longitudinal  Response 


Balloon  Displa  'ements 
From  Equilibri'.m 

Cable  Tension 
At  Equilibrium 

Max  Cable 
Tension 

Case 

Winch 

Alt. 

Range 
Down/ Up 

Alt. 

Pitch 

At 

Winch 

At 

Bridle 

At 

Winch 

At 

Bridle 

No. 

Ft. 

Ft. 

Ft. 

Deg. 

Lb. 

Lb. 

Lb. 

Lb., 

3 

0.  0 

-52  /73 

82 

2.9 

1650 

2210 

2650 

3210 

17 

5000 

-20  /47 

24 

3.  1 

1940 

2210 

3670 

3940 

(10)  Effects  of  Different  Operational  Altitudes.  The  effects  of  three 
different  operational  altitudes  on  tho  response  of  BJ  type  balloons  are  given 
in  Table  XXII.  The  design  operational  altitudes  and  corresponding  balloon 
volumes  are  5,000  ft.  and  46,000  cu.  ft.,  10,000  ft.  and  60,000  cu.  ft.,  and 
20,000  ft.  and  500,000  cu.  ft.  The  balloons  are  in  equilibrium  with  steady 
head  winds  which  vary  according  to  wind  velocity  altitude  profile.  Longitudinal 
ten  second  duration  gusts  of  15.7  ft/  sec  ,  20.2  ft/sec,  and  27.3  ft/sec  are 
applied  respectively  to  the  5,000,  10,000,  and  20,000  ft.  operational  altitude 
configurations  in  the  longitudinal  cases.  Lateral  three  second  ramp  rise  time 
steady  gusts  of  9  ft/sec,  12  ft/sec,  and  16  ft/sec  are  applied  respectively  to 
the  5,000,  10,000,  and  20,000  ft.  operational  altitude  configurations  in  the 
lateral  cases.  Table  XXI  reveals  that 

(1)  the  range  displacements  from  equilibrium  increase  with  an 

increase  in  operational  altitude 

(2)  pitch  displacements  from  equilibrium  decreases  with  an 
increase  in  operational  altitude 

(3)  the  lateral  responses  show  little  change  due  to  operational 
altitude 

(4)  the  tether  load  increases  wirh  altitude  as  expected. 


Table  XXI.  Effect  of  Various  Design  Operational  Altitudes 


From  Equilibrium 


Casv;  Volume 


8  46000 

3  60000 

9  500000 


Gust* 

Range 

Vel. 

Down/Up 

Alt. 

Ft/Se< 

:  Ft. 

Ft. 

15.  7 

-12  /38 

17 

20.  2 

-52  /73 

82 

27.  3 
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Cable  Tension 
At  Equilibrium 


Max  Cable 
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B.  Balloon  Structural  Design  Performance  Parameters 


The  balloon  structural  design  performance  parameters  of  interest  are 
the  tensions  at  the  winch  and  bridle  and  the  translational  and  rotational 
accelerations  of  the  balloon  center  of  gravity.  Tables  XXIV  and  XXV  list  the 
maximum  variation  of  these  parameters  from  equilibrium  for  selected  longitudinal 
and  lateral  cases.  From  a  design  limit  load  point  of  view,  the  values  found 
in  Table  XXII  represent  the  structural  design  parameters  for  the  three  balloon 
types  subjected  to  longitudinal  gusts  of  20.2  ft/sec  and  10  second  duration  and 
lateral  steady  gusts  of  12  ft/sec. 


Table  XXII.  Structural  Design  Parameters 
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.4 

.  7 

.  5 
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C.  Balloon  Instrumentation  Design  Performance  Parameters 


The  balloon  instrumentation  design  performance  parameters  of  interest 
in  designing  an  "on-board"  instrumentation  package  are  as  follows: 

(1)  displacements  (horizontal,  vertical,  and  lateral) 

(2)  velocities  (horizontal,  vertical,  and  lateral) 

(3)  angular  displacements  (pitch,  yaw,  and  roll) 

(4)  angular  velocities  (pitch,  yaw,  and  roll) 

Tables  XXIV  and  XXV  list  the  maximum  variation  of  these  parameters  from 
equilibrium  for  selected  longitudinal  and  lateral  cases. 

From  an  instrumentation  design  point  of  view,  the  values  found  in 
Table  XXIII  represent  the  instrumentation  design  parameters  for  the  three 
balloon  types  subjected  to  longitudinal  gusts  of  20.2  ft/sec  and  10  second 
duration  and  lateral  steady  gusts  of  12  ft/sec. 


Table  XXIII.  Instrumentation  Design  Parameters 
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COMPARISON  OF  BALLOON  TYPES 


The  BJ  balloon,  with  ram  air  filled  ballonet  and  fins,  provides  the 
smallest  size  tethered  balloon  system  to  fly  a  lOOo-pound  payload  at 
10,000  feet.  However,  it  has  the  greatest  downrange  displacement  from 
the  ground  tether  point.  Equilibrium  tether  tensions  are  the  lowest. 
Longitudinal  excursions  around  equilibrium  are  moderate,  but  lateral  dis¬ 
placements  due  to  side  gusts  are  substantial  as  compared  to  the  Vee-type 
balloon. 


Generally  the  Vee  balloon  investigated  has  a  somewhat  greater 
angle-of-attack  excursion  about  static  equilibrium  for  a  given  input,  but 
with  the  high  aerodynamic  lift  configuration,  results  in  less  overall  downrange 
displacement  from  the  tether  point.  Lateral  excursion  is  a  minimum  for  this 
balloon  type.  Overali,  the  Vee  offers  the  least  excursion  with  wind  variation 
but  at  the  expense  of  relatively  high  tether  tension.  Elastic  devices  can  be 
built  into  the  suspension  system  to  limit  angle  of  attack  and  tether  tension. 
However,  analysis  of  elastic  suspension  is  beyond  the  scope  of  the  reported 
analysis.  The  Vee-type  balloon  then  offers  a  configuration  for  applications 
where  the  tethered  system  should  be  nearer  vertical  with  respect  to  the  tether 
point,  and  where  displacement  from  the  tether  point  should  be  minimized. 

The  GAG  single-hull  balloon  provides  a  greater  aerodynamic  lift  and 
aerodynamic  lift-to-drag  ratio  and  less  downrange  displacement  from  the 
tether  point  than  the  EJ  balloon,  ind  not  much  greater  than  that  for  the  Vee 
for  the  10,000-foot  altitude  system  (4950  feet  as  compared  to  7377  feet  for 
the  BJ  ana  4594  feet  for  the  Vee).  The  equilibrium  cable  tension  is  moderate 
(3800  pounds).  Longitudinal  motions  of  this  balloon  are  substantially  less  for 
given  input  wind  disturbances  than  for  the  other  types.;  Lateral  displacements 
are  comparable  to  those  for  the  BJ  balloon  for  an  equal  disturbance.  The  first 
longitudinal  and  lateral  modes  of  GAC  single-hull  ballot>n  type  are  notably  more 
stable  (set  Figures  31  and  51  of  Reference  2)  than  the  other  balloon  types 


Structural  and  onboard  payload  design  parameters  for  particular  wind 
gust  inputs  are  listed  in  Tables  XXIV  and  XXV. 

3.  INFLUENCE  OF  TRIM  ANGLE  ON  LONGITUDINAL  AND  LATERAL  DYNAMIC  BEHAVIOR 

Considering  the  BJ  balloon,  longitudinal  dynamic  response  dara 
indicates  that  trim  angle-of-attack  does  not  significantly  change  the  response 
of  the  balloon  in  longitudinal  translation  and  rotation.- 

In  view  of  the  lateral  dynamic  response,  analysis  indicates 
that  the  two  lower  modes  of  motion  are  excited  (i.e,,  pendulum  motion  of 
the  balloon  and  tether  as  a  whole,  and  balloon  roll),.  Trim  angle  change  does 
not  substantially  influence  these  modes  for  the  BJ  balloon.  However,  the 
vertical  location  of  the  confluence  point  affects  the  second  mode  (balloon  roll) 
The  apex  closer  to  the  balloon  results  in  greater  damping  of  this  mode  while 
not  greatly  affecting  the  dat.ping  of  the  first  mode  (Figure  36  of  Reference  2), 
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Longitudinal  dynamic  response  of  the  Vee  balloon  is  somewhat  more 
sensitive  to  design  trim  angle-of-attack.  A  tethered  balloon  system, 
employing  the  Vee  balloon,  should  be  designed  to  fly  at  moderate  angle-cf- 
attack  (say  7  degrees)  inasmuch  as  longitudinal  instability  tendencies  °vist 
at  low  angle-of-attack.  This  tendency  is  indicatec  in  the  linearized  stability 
analysis  by  the  lack  of  damping  in  the  roots  of  the  characteristic  equation  for 
the  lowest  frequency  mode  of  motion.  From  the  stability  analysis,  first  and 
second  longitudinal  and  lateral  modes  of  motion  for  the  GAC  single  hull  balloon 
are  stable  and  well  damped  (Figures  31  and  51  of  Reference  2).  The  first  two 
longitudinal  modes  are,  respectively,  a  pitching  motion  of  tether  and  balloon, 
and  a  pitching  of  the  balloon.  The  first  two  lateral  modes  are  a  coupled  balloon 
yaw  and  lateral  link  rotation,  and  a  coupled  balloon  yaw  and  roll  and  lateral 
link  rotation.  These  modes  of  motion  are  relatively  insensitive  to  trim  angle 
of  attack  and  vertical  location  of  the  bridle  confluence  point.  Dynamic  response 
data  also  verify  that  longitudinal  linear  and  angular  displacements  from 
equilibrium  are  relatively  insensitive  to  trim  angle-of-attack  changes, 

4.  EFFECT  OF  TAIL  SIZE  ON  LONGITUDINAL  AND  LATERAL  DYNAMIC  BEHAVIOR 

Tail  size  for  the  BJ  balloon  was  varied  by  proportionally  changing 
all  dimensions  to  obtain  tail  area  of  81%,  100%  and  144%  of  the  nominal 
design  values.  Stability  analysis  indicates  that  tail  size  increase  to  144%  is 
required  to  give  longitudinal  stability  for  the  first  oscillatory  mode, but  that 
all  tail  sizes  provide  lateral  stability  in  the  first  lateral  mode.  Also,  the 
tether  dynamic  load  factor  reduces  with  larger  tail  sizes  because  of  reduced 
angle-of-attack  excursions.  However,  increased  tail  sizes  increase  the  side 
displacements  somewhat.  It  appears  that  100%  tails  should  be  adequate  since 
near  neutral  stability  is  present  about  the  equilibrium  point  for  the  first 
longitudinal  mode,  and  dynamic  response  data  does  not  show  displacement  diver¬ 
gence  for  any  of  the  tail  sizes  for  the  180  seconds  of  response  that  were 
calculated. 

The  Vee  balloon  investigation  was  confined  to  dete’-mining  the  effect 
of  increasing  vertical  tail  area  only.  Consequently,  longitudinal  dynamic 
characteristics  are  not  greatly  affected.  The  vertical  tails  on  the  lower  side 
of  the  two  hulls  are  increased  in  area  from  the  100%  nominal  size  to  300%. 

An  examination  of  the  roots  of  the  characteristic  equations  for  lateral  stability 
(Figure  49  of  Reference  2)  indicates  that  the  larger  tails  give  somewhat 
greater  damping  and  lower  frequency  for  the  first  mode,  and  an  inverse  effect 
for  the  second  mode.  The  lateral  dyi.amic  response  data  indicate  that  lateral 
displacements  of  yaw,  roll,  and  side  motion  are  all  reduced  with  the  larger 
tail  sizes.  Tail  areas  of  200%,  of  nominal  then  appears  to  be  most  desirable. 

5.  EFFECT  OF  ALTITUDES  INTERMEDIATE  TO  TEE  DESIGN  ALTITUDE  ON  BEHAVIOR 
OF  THE  BJ  TETHERED  BALLOON  SYSTEM 

As  the  3J  tethered  balloon  syr'.em  designed  for  10,000  feet  is  brought 
to  lower  and  lower  altitudes,  the  wind  speed  it  sees  is  reduced  and  the  trim 
angle  of  attack  changes  inasmuch  as  the  suspension  bridle  geometry  remains 
unchanged.  Generally,  the  longitudinal  modes  of  motion  become  less  damped 
at  the  lower  altitudes.  Dynamic  response  data  indicates  that  longitudinal 
response  is  diminished,  as  might  be  expected,  with  the  lower  steady  winds  and 
gusts  used  in  the  analysis. 
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The  lateral  stability  analysis  reveals  an  unstable  lateral  first  mode 
consisting  of  coupled  balloon  yaw  and  lateral  displacement  of  the  1000-foot 
altitude.  Generally,  damping  of  all  lateral  modes  is  reduced  substantially 
at  the  lower  altitudes,  A  divergent  lateral  motion  is  revealed  in  the  dynamic 
response  analysis  at  1000  feet  altitude  for  the  10,000-foot  design.  Comparing 
dynamic  response  and  stability  data,  it  is  concluded  that  the  unstable  first 
lateral  mode  consists  of  a  coupled  yaw  lateral  displacement.  This  mode 
oscillates  at  a  frequency  which  is  approximately  equal  to  an  inverted  pendulum 
mode  where  the  frequency  is  established  by  the  pendulum  length  and  the 
restoring  forces  consisting  of  the  net  buoyancy  and  the  aerodynamic  forces. 

The  divergent  motion  was  not  revealed  at  the  4000-foot  altitude.  Design 
features  which  improve  the  damping  of  this  mode,  such  as  larger  tails  and/or 
automatically  actuated  stabilizing  aerodynamic  surfaces,  are  required.  A  more 
detailed  investigation  and  analysis  is  required  to  make  specific  recommendations 
to  improve  the  dynamic  characteristics  exposed  by  this  exploratory  analysis. 

6.  EFFECTS  OF  REDUCED  WIND  ON  THE  BEHAVIOR  OF  THE  BJ  TETHERED  BALLOON  SYSTEM 

The  longitudinal  and  lateral  dynamic  characteristics  consisting  of 
frequencies  and  damping  qualities  change  considerably  as  the  wind  is  reduced.  The 
trim  angle  of  attack  steadily  decreases  as  wind  .decreases.  Generally,  damping 
of  the  tethered  balloon  system  also  reduces  substantially  as  wind  is  reduced. 

As  steady  winds  and  wind  gusts  reduce  yaw  response  is  unchanged  but  other 
balloon  motion  responses  reduce  as  might  be  expected  . 

7.  NOLARO  VS  AMGAL  TETHERS 

The  BJ  balioon  system  employing  a  NOLARO  type  tether  offers  the 
advantage  of  a  smaller  balloon  to  meet  the  requirements  of  supporting  a  1000- 
pcund  payload  at  10,000  feet.  Dynamic  response  of  the  balloon  longitudinal 
and  lateral  degrees  of  freedom  are  not  greatly  influenced  by  the  different  tethers 
with  the  exception  of  lateral  displacements.  Greater  lateral  displacements 
result  with  the  system  employing  AMGAL.  This  may  be  attributed  to  the  lesser 
lateral  damping  of  the  smaller  diameter  AMGAL  tether.  For  20,000  foot 
altitude  designs,  a  single  balloon  tethered  system  can  not  be  achieved  with  a 
reasonable  size  balloon  using  AMGAL  tether  material.  Considering  these 
factors,  the  NOIARO  type  tether  generally  will  provide  a  more  favorable 
tethered  balloon  system. 

8.  EFFECT  OF  PAYLOAD  LOCATION 

The  configuration  with  payload  located  at  the  confluence  point  of  the 
bridle  has  somewhat  greater  longitudinal  response  than  that  with  the  payload 
located  at  the  underside  of  the  balloon.  However,  longitudinal  response 
changes  are  relatively  minor  and  need  not  influence  i  decision  on  payload 
location.  Lateral  dynami"  response  characteristics  were  not  determined  for 
this  design  variation. 

9.  EFFECT  OF  WINCH  ALTITUDE  LOCATION 

In  this  comparison  the  BJ  balloon  is  still  flown  at  10,000  feet  altitude, 
M.S.L.,  but  the  winch  is  at  5000  feet  M.S.L.  Dynamic  response  is  somewhat 
larger  in  the  longitudinal  plane  with  the  winch  at  5000  feet.  Lateral  dynamic 
response  characteristics  were  not  determined  for  this  design  variation.  Lateral 
stability  analysis  indicates  that  one  of  the  overdamped  modes  becomes  oscillatory 
when  the  winch  is  located  at  5000  feet. 
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10.  BJ  TETHERED  BALLOON  SYSTEMS  DESIGNED  FOR  VARIOUS  ALTITUDES 


Tethered  balloon  system  designs  for  altitudes  of  5000,  10,000  and 
20,000  feet  have  been  defined.  As  commented  elsewhere,  it  was  not  possible 
to  obtain  a  reasonably  sized  single  balloon  system  with  AMGAL  tether  for  the 
20,000  foot  condition.  Dynamic  response  analysis  indicates  that  longitudinal 
displacements  are  somewhat  reduced  for  higher  altitude  designs  and  lateral 
displacements  are  somewhat  increased.  However,  substantial  changes  have  not  been 
observed.  Dynamic  load  factors  for  tether  cable  design  are  significantly  less 
for  the  higher  altitude  design. 

Longitudinal  stability  analysis  indicates  that  the  second  and  third 
oscillatory  modes,  which  are  balloon  pitching  modes,  reduce  in  frequency  for 
the  higher  altitude  designs.  The  lower  frequencies  might  be  somewhat 
attributed  to  the  larger  pitching  inertia  associated  with  the  larger  balloons. 

It  is  recommended  that  the  mathematical  techniques  developed  be 
investigated  further  by  comparing  analytical  predictions  of  the  tools  with 
experimental  data.  After  this  correlation  has  been  completed,  more  compre¬ 
hensive  analysis  of  tethered  balloon  systems  of  specific  interest  with  wind 
conditions  and  gusts  to  be  expected  in  operations  can  be  conducted. 
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APPENDIX  A 

EQUATIONS  OF  MOTION 
FOR  DYNAMIC  SIMULATION 
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SECTION  I 
INTRODUCTION 


In  Appendix  A  of  the  Second  Scientific  Report*,  the  equations 
of  motion  of  a  tethered  balloon  were  derived  in  three  dimensions. 
The  initial  derivation  assumed  that  angular  velocities  were  small 
so  that  products  of  angular  velocities  could  be  neglected.  This 
assumption  reduced  the  complexity  of  the  equations  (especially 
the  lateral  equations)  by  a  considerable  amount,  but  each 
equation  remained  nonlinear  and  coupled  in  each  degree  of  freedom. 
If  is  assumed  that  the  reader  has  read  Ref.  2  Appendix  A,  Section 
II  and  has  a  clear  understanding  of  the  derivation  of  the  equation 
of  motion,  for  this  discussion  will  begin  with  the  results  of  that 
section. 

Many  of  the  terms  used  in  this  derivation  are  defined  in  the 
listing  which  follows. 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

AA(4,  4) 

% 

A  four  by  four  array  of  variables 
used  to  define  incremental 
velocities  in  the  numerical 
integration 

rad/sec 

AAG (8) 

An  array  of  eight  variables 
signifying  the  angle  of  the 
gust  to  the  horizon  corre¬ 
sponding  to  TTG (8) 

deg 

AALP (16) 

An  array  of  sixteen  variables 
signifying  the  angle-of-attack 
of  the  balloon 

rad 

AALPD (16) 

An  array  of  sixteen  variables 
signifying  the  angle-of-attack 
of  the  balloon 

deg 

AALT (8 ) 

An  array  of  eight  variables 
signifying  altitudes  of  steady 
state  wind  velocities 

ft 

AG 

ag 

Gust  angle  off  the  horizon 
interpolated  from  AAG (8)  and 

TTG { 8 ) 

deg 

AGR 

ag 

Gust  angle  off  the  horizon 

rad 

ALPB 

“b 

Angle-of-attack  of  balloon 

rad 

ALPBD 

& 

B 

d  afl/dt 

rad/sec 

ALPBDD 

“b 

d  aB/dt 

deg/sec 

ALPBDE 

a 

B 

Angle-of-attack  of  balloon 

deg 

ALPSL 

A  ratio  of  two  angle-of-attack 
differences  used  in  inter¬ 
polation  of  aerodynamic 
coefficients  of  the  balloon  at 
some  a 

ALTSL 

A  ratio  of  two  altitude  differ¬ 
ences  used  in  interpolation  of 
steady  state  wind  velocities 

AXB 

^XB 

The  inertial  force  acting  on  the 
balloon  along  the  lateral  axis 
of  the  balloon 

slug-f t/ 
2 

sec 

AYB 

Uyb 

The  inertial  force  acting  on 
the  balloon  along  the  longitudinal 
axis  of  the  balloon 

slug-f t/ 
sec^ 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

AZB 

^ZB 

The  inertial  force  acting  on 
the  balloon  along  the 
vertical  axis  of  the  balloon 

slug-ft/ 

sec^ 

BBET ( 8 ) 

An  array  of  eight  variables 
signifying  the  sideslip 
angle  of  the  balloon 

rad 

BBETD ( 8 ) 

An  array  of  eight  variables 
signifying  the  sideslip 
angle  of  the  balloon 

deg 

BETJ 

0 

Sideslip  angle  of  balloon 
(positive-wind  to  the  right 
of  balloon's  centerline) 

rad 

BETBD 

• 

d  0/dt 

rad/sec 

BETBDD 

• 

6 

d  0/dt 

deg/sec 

BETBDE 

6 

Sideslip  angle  of  balloon 
(positive-wind  to  the  right 
of  balloon's  centerline) 

deg 

BETSL 

A  ratio  of  two  sideslip  angle 
differences  useu  in  inter¬ 
polation  of  aerodynamic 
coefficients  of  the  balloon 
at  some  0 

C  (3) 

C 

r 

Non-dimensional  center-of- 
pressure  of  "r"  th  link 

CCDB (16) 

An  array  of  sixteen  variables 
signifying  drag  coefficients 
of  the  balloon  corresponding 
to  AALPD (16) 

CCLB  (16) 

An  array  of  sixteen  variables 
signifying  lift  coefficients 
of  the  balloon  in  the  longi¬ 
tudinal  program  corresponding 
to  AALPD (16) 

CCLB  ( 8 ) 

An  array  of  eight  variables 
signifying  roll  moment  coef¬ 
ficients  of  the  balloon  in  the 
lateral  proqram  corresponding 
to  BBETD (8)  ' 

CCLL3 (16) 

An  array  of  sixteen  variables 
signifying  lift  coefficients 
of  the  balloon  in  the  lateral 
program  corresponding  to  AALPD (16) 

A-5 


4 


FORTRAN  STANDARD 
CCMB (16) 


CCNB (8) 


DESCRIPTION  UNITS 


An  array  of  sixteen  variables 
signifying  pitch  moment 
coefficients  of  the  balloon 
corresponding  to  AALPD(16) 

An  array  of  eight  variables 
signifying  yaw  moment 
coefficients  of  the  balloon 
corresponding  to  BBETD ( 8 ) 


CDB 

CDB 

CD  C 

CDC 

CDTDB 

CD0E 

CGAMB 

C  ^B 

CLB 

CLB 

CLB 

C  S.B 

CLLB 

CLB 

CLPHDB 

c  • 

£<J)B 

CLPSDB 

C4.i|jB 

CLTDB 

CL§B 

CMC 

CmB 

CMTDB 

CNB 

CnB 

CNPHDB 

Cn$9 

Drag  coefficients  of  balloon 

Drag  coefficient  of  link 
(infinite  cylinder) 

Balloon  drag  coefficient  due 
to  pitch  velocity  rad 

cos  yb 

Lift  coefficient  of  balloon 
in  longitudinal  program 

Balloon  roll  moment  coefficient 
in  lateral  program 

Lift  coefficient  of  balloon  in 
lateral  program 

Balloon  roll  moment  coefficient 
due  to  roll  velocity  rad 

Balloon  roll  moment  coefficient 
due  to  yaw  velocity  rad 

Balloon  lift  coefficient  due 
to  pitch  velocity  rad 

Balloon  pitch  moment 
coefficient 

Balloon  pitch  moment  coef¬ 
ficient  due  to  pitch  velocity  rad 

Balloon  yaw  moment  coefficient 

Balloon  yaw  moment  coefficient 
due  to  roll  velocity  rad 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

CNPSDB 

c  • 
ni^B 

Balloon  yaw  moment  coef¬ 
ficient  due  to  yaw 
velocity 

rad  d 

COM1 

Input  data  defining  com¬ 
puter  run 

CPHI 

C<f> 

COS  ()> 

CPHI2 

C2(J) 

2  x 

cos  <j> 

CPSI 

C  ij> 

cos  ip 

CSIG  (3) 

coi 

cos  =  1,  2,  3) 

CSIG2 (3) 

C2a . 
l 

cos2oi  x(i  =  1,  2,  3) 

CTHE 

C0 

cos  0 

CTHE2 

c2e 

cos2  0 

CTPX ( 3 ) 

C(0+5£) 

cos  (0  +  q)  ,  (i  =  1,  2,  3) 

CTPX2 (3) 

c2(e+q> 

cos2  (0  +  £;i)  t  (i  =  1,  2,  3) 

CXI (3) 

cos  ,(i  =  1,  2,  3) 

CXI2 (3) 

cos2  £i  ,  (i  =  1, 2,3) 

CXMX ( 3 , 3) 

cos(Ci  -Cj)  ,(i  =  1,2,3,;  j=  1, 

2,3) 

CYB 

CYB 

Side  force  coefficient 

CYPHDB 

C**B 

Side  force  coefficient  of 
balloon  due  to  roll  velocity 

:ad_1 

CYPSDB 

Side  force  coefficient  of 
balloon  due  to  yaw  velocity 

rad 

DB 

dB 

Aerodynamic  reference  length 
of  balloon  (¥■*•/' ) 

ft 

DC 

Diameter  of  tether 

ft 

jCH(3) 

DCHi 

Horizontal  aerodynamic  force 
acting  on  the  center-of- 
pressure  of  the  "i"th  link 
(i  =  1,2,3) 

lbf 

A-7 
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FORTRAN 

DCS (3) 

DCV ( 3 ) 

DD (4,4) 

DRAGB 

DT 

DTP 

DTP1 

DTP  2 

DTI 

DT2 

DYPRB 

DYPRH (4 ) 


STANDARD  DESCRIPTION  UNITS 


CSi 


D 


CVi 


DRAGB 


^B 


q 


Hi 


Lateral  aerodynamic  force 
acting  on  the  center-of- 
pressure  or  the  "i"th  link, 

(i  =  1/2,3) 

lbf 

Vertical  aerodynamic  force 
acting  on  the  center-of- 
pressure  of  the  "i"th  link/ 

(i  =  1,2,3) 

lbf 

A  four  by  four  array  of 
variables  signifying  the 
coefficients  of  the  second 
derivatives  in  the  equations 
of  motion 

slug-ft 

Aerodynamic  drag  acting  on 
balloon 

Integration  time  increment 

sec 

Number  of  integrations  be¬ 
tween  data  printouts 

Number  of  integrations  be¬ 
tween  data  printouts  when 

DT  =  D'.?l 

Number  of  integrations  be¬ 
tween  data  printouts  when 

DT  =  DT2 

Integration  time  increment 
when  T  <  TDTC 

sec 

Integration  t^me  increment 
when  T  >,  TDTC 

sec 

Effective  dynamic  pressure 
acting  on  c.g.  of  balloon 

lbf/ft2 

Effective  dynamic  pressure 
acting  on  the  "i"th  link 
hinge  ,  (i  =  1,2, 3, 4).  The 
first  hinge  is  at  the  winch 

lbf/tt* 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

EE  (4) 

An  array  of  four  variables 
signifying  the  generalized 
forces  plus  terms  containing 
products  of  angular  velocities 
if  this  opt ‘on  is  chosen. 

See  FF 

ft-lb 

EEPS (3,3) 

e.  . 

ij 

A  three  by  three  array  re¬ 
presenting  the  transformation 
matrix  between  the  balloon 
body  unit  vectors  and  the 
inertia  frame  unit  vectors 

FBHA 

F 

BHA 

Horizontal  aerodynamic  forces 
acting  on  balloon 

lbf 

FBSA 

F 

BSA 

Lateral  aerodynamic  forces 
acting  on  balloon 

lbf 

FBVA 

F 

BVA 

Vertical  aerodynamic  forces 
acting  on  balloon 

lbf 

FCH  (3) 

F 

CHi 

Horizontal  component  of 
tension  acting  on  the  ''i"th 
link  ,(i  =  1,2,3) 

lbf 

FCN (3) 

F 

CNi 

Aerodynamic  force  acting 
normal  to  " i"th  link, 

(i  =  1,2,3) 

lb* 

FCS  (3) 

F 

CSi 

Lateral  component  of  tension 
acting  on  the  "i"th  link, 

(i  =  1,2,3) 

lb£ 

FCV ( 3 ) 

F 

CVi 

Vertical  component  of  tension 

acting  on  the  "i"t.h  link 

(i  -  1,2,3)  lbf 


FF  Input  variable  which  gjves 

the  option  of  including 
inertia  terms  containing  pro¬ 
ducts  of  angular  velocities 
in  the  longitudinal  program. 
If  these  terms  are  not  wanted 
set  FF  =  0.  .  If  the  terms  are 
to  be  included,  set  FF  to  any 
number  except  "0.". 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

FPHI 

F, 

<P 

Generalized  force  acting  on 
balloon's  roll  degree-of- 
f reedom 

ft-lbf 

FPSI 

F  , 

if 

Generalized  force  acting  on 
balloon's  yaw  degree-of- 
freedom 

f t-lbf 

FSIG (3) 

F*. 

ai 

Generalized  force  acting  on 
the  "i"th  link's  yaw  degree- 
of-freedom; (i  =  1,2,3) 

f t-lbf 

FTHE 

pe 

Generalized  force  acting  on 
balloon's  pitch  degree-of- 
freedom 

ft-lbf 

FWH 

Horizontal  component  of  tether 
tension  acting  at  winch 

lbf 

FWS 

Lateral  component  of  tether 
tension  acting  r.t  winch 

lbf 

FWV 

Vertical  component  of  tether 
tension  acting  at  winch 

lbf 

FXI  (3) 

F5i 

Generalized  force  acting  on 
the  "i"th  link's  pitch  degree- 
of-freedomy (i  -  1,2,3) 

ft-lb^ 

G 

g 

Acceleration^of  gravity  = 

32.17  ft/sec 

ft/sec 

GAMB 

Y3'  rB 

I 

In  the  longitudinal  program 

Yb  is  the  relative  flight  path 
angle  of  the  balloon  in  the 
longitudinal  plane;  in  the 
lateral  program  Fp  is  the 
relative  flight  path  angle  of 
the  balloon  in  the  horizontal 
plane 

rad 

GAMBD 

y  r 

'b,  b 

d  Yg/dt  or  a  Fg/dt 

rad/sec 

GAMBDE 

y  r 

B,  !B 

See  GAMB 

deg 

GAi'lBDD 

.  « 

y  r 

B,  B 

dYB/dt  or  d  F^dt 

deg/sec 

IXB 

IXB 

Apparent  pitch  moment  of 
inertia  of  balloon 

slug-f t 

A  - 1 0 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

IYB 

XYB 

Apparent  roll  moment  of 
inertia  of  balloon 

slug-f t 

IYZB 

IYZB 

Apparent  product  of.  inertia 
of  balloon  with  respect  to 
roll  and  yaw  body  axis 

slag-ft 

IZB 

IZB 

7  'parent  yaw  moment  of 
i.  ertia  of  balloon 

slug-ft 

L 

£ 

Length  of  one  tether  link 

ft 

LIFTB 

LIFTB 

Aerodynamic  lift  acting  on 
balloon 

ibf 

LS 

Ls 

Buoyant  lift  of  balloon  = 
weight  of  displaced  air 

lbf 

M 

m 

Mass  of  one  tether  link 

slugs 

MAL 

Added  mass  of  balloon  along 
its  longitudinal  axis 

slugs 

MAS 

Added  mass  of  balloon  along 
its  lateral  axis 

slugs 

MAV 

Added  mass  of  bal'oon  along 
its  vertical  axis 

slugs 

ML 

ml 

Apparent  mass  of  balloon 
along  its  longitudinal  axis. 
Apparent  mass  equals  total 
static  mass  associated  with 
the  balloon  plus  air  mass  being 
accelerated  in  a  particular 
direction. 

[ 

slugs 

MPHIB 

M4>b 

Aerodynamic  rolling  moment 
acting  on  balloon 

f t-lbf 

MPL 

mPL 

Mass  of  payload-located  at 
bridle  tether  connection 

slugs 

MPSIB 

M’^B 

Aerodynamic  yawing  moment 
acting  on  balloon 

ft-lbf 

MS 

r,s 

Apparent  mass  of  balloon 
along  its  lateral  axis 

slugs 

MTHEB 

MeB 

Aerodynamic  pitching  moment 
acting  on  balloon 

ft-lb, 

f 

A  - 1 1 
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FORTRAN 

STANDARD 

DESCRIPTION 

MV 

“v 

Apparent  mass  of  balloon 
along  its  vertical  axis 

PHI 

<J> 

Roll  angle  of  balloon 
(positive-right  wing  down) 

PHID 

$ 

d<j>/dt 

PHIDD 

V 

2  2 
d  ?/dtz 

PHIDDD 

2  2 
dz<))/dt 

PHIDDE 

d<f>/dt 

PT1IDEG 

Roll  angle  of  balloon 
(positive-right  wing  down) 

PHID2 

4>2 

PR 

Atmospheric  pressure  (not 
used  in  program) 

PSI 

¥ 

Yaw  angle  of  balloon 
(positive-nose  to  right) 

PSID 

* 

dijj/dt 

PSIDD 

¥ 

2  2 
d  >/dt 

PSIDDD 

•t 

2  2 
d>/dt 

PSIDDE 

* 

dtp/dt 

PSIDEG 

'P 

Yaw  angle  of  balloon 
(positive-nose  to  right) 

PSID2 

*2 

^2 

RHOB 

PB 

Air  density  at  Z 

O 

RHOC  ( 3 ) 

pci 

Air  density  at  Z  .  ,(i=l,2,3) 

RHOCH ( 3 ) 

pHi 

Air  density  at  Zc^^((i=l ,  2 , 3 , 4 

RJA 

R.. 

da 

Distance  along  longitudinal 
axis  of  balloon  from  aero¬ 
dynamic  reference  center  to 

bridle  apex  (positive 
toward  nose) 


UNITS 

slugs 

rad 

rad/sec 

rad/sec2 

deg/sor2 

deg/sec 

deg 

2  2 

rad  /sec 

lb/ft2 

rad 

rad/sec 
rad/sec2 
deg/ sec 2 
deg/sec 

dc: 

2  2 

rad  /sec 

slug/ft2 

slug/ft2 

)  slug/ft2 

ft 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

RJB 

RjB 

Distance  along  longitudinal 
axis  of  balloon  from  center 
of  buoyancy  to  bridle  apex 
(positive  toward  nose) 

ft 

RJG 

Rjg 

Distance  along  longitudinal 
axis  of  balloon  from  center 
of  gravity  to  bridle  apex 
(positive  toward  nose) 

ft 

RJM 

Rjm 

Distance  along  longitudinal 
axis  of  balloon  from 
dynamic  mass  center  to 
bridle  apex  (positive 
toward  nose) 

ft 

RKA 

RkA 

Distance  along  vertical 
axis  of  balloon  from  aero¬ 
dynamic  reference  center  to 
bridle  apex  (positive  up) 

ft 

RKB 

RkB 

Distance  along  vertical  axis 
of  balloon  from  center  of 
buoyancy  to  bridle  apex 
(positive  up) 

ft 

RKG 

K 

kg 

Distance  along  vertical 
axis  of  balloon  from  center 
of  gravity  to  bridle  apex 
(positive  up) 

ft 

RKM 

Rkm 

Distance  along  vertical  axis 
of  balloon  from  dynamic  mass 
center  to  bridle  apex 
(positive  up) 

ft 

SB 

SB 

Aerodynamic  reference  area 
of  balloon  (Vt) 

ft2 

SC  { 3) 

Sci 

Aerodynamic  reference  area 
of  one  tether  link  =  (l)  (dc) 

(i  =  1/2,3)  ' 

n2 

SGAMB 

S^B 

sin  yb 

SIG  { 3) 

a . 

i 

Yaw  angle  of  "i"th  link; 

(i  =  1,2,3)  (positive-counter¬ 
clockwise  looking  down  range) 

rad 

SIGD ( 3) 

*i 

da i/dt  t  (i  =  1,2,3) 

rad/sec 

A  - 1  3 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

SIGDD(3) 

•  % 

cr. 

i 

dV/dt2  ,(i  =  1,2,3) 

rad/sec2 

SIGDDD{3) 

«  . 
o. 
i 

d2ai/dt2  ;(i  =  1,2,3) 

deg/sec2 

SIGDDE (3) 

• 

0  . 

1 

da^/dt  ,{i  -  1,2,3) 

deg/sec 

SIGDEG (3) 

0  . 

1 

Yaw  angle  of  "i"th  link/ 

(i  =  1,2,3)  (positive- 
counterclockwise  looking 
down  range) 

deg 

SIGD2 (3) 

d.2 

1 

V  ,(i=  1,2,3) 

rad2/sec" 

SPHI 

s<p 

sin  4> 

SPHI2 

S 2  4> 

sinz4> 

SPSI 

sin  ijj 

SSIG (3) 

Sa . 
i 

sin  i (i  =  1,2,3) 

STHE 

S0 

sin  0 

STHE2 

S20 

sin2  0 

STPX (3) 

s ( 0  +ci) 

sin  (0+E^)  ,(i  =  1,2,3) 

STPX2 (3) 

S2(Q+^±) 

sin  ^0+CJ  ;  (i  =  1,2,3) 

SXI (3) 

sin  ,  (i  =  1,2,3) 

SXI2 (3) 

s\ 

sin2Cf  ,  (i  =  1,2,3) 

SXMX (3,3) 

s(5i-5j) 

sintSj-Sj)  /  (i=l, 2, 3;  j=l, 2, 3) 

S2THE 

S20 

sin  20 

T 

TDTC 

t 

Flight  time 

Flight  time  at  which  DT  is 
changed  from  DTI  to  DT2  and 

DTP  is  changed  from  DTPl  to 

DTP  2 

sec 

sec 

TENS (3) 

T  . 

l 

Tension  in  tether  at  top  of 
" i " th  link  ;(i  =  1,2,3) 

lbf 

TENSW 

T 

w 

Tension  in  tether  at  winch 

lbf 

A  - 14 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

TETH 

Total  length  of  tether 

ft 

THE 

0 

Pitch  angle  of  balloon 
(positive-nose  up) 

rad 

THED 

• 

0 

d0/dt 

rad/sec 

THEDD 

•  C 

0 

d20/dt 

rad/sec2 

THEDDD 

•  • 

0 

d2  0/dt 

deg/sec2 

THEDDE 

• 

0 

d0/dt 

deg/sec 

THEDEG 

0 

Pitch  angle  of  balloon 
(positive-nose  up) 

deg 

THED2 

02 

02 

rad2/sec 

THEODE 

0 

o 

Equilibrium  pitch  angle 
of  balloon 

deg 

TSL 

A  ratio  of  two  time 
differences  used  in  inter¬ 
polation  of  gust  velocities 
and  gust  angles 

TTG (8) 

An  array  of  eight  variables 
signifying  the  time  history 
of  the  gust 

sec 

TTT 

Flight  time  at  which 
computer  simulation  ends 

sec 

VBR 

VBR 

Relative  velocity  of 
balloon's  c.g.  with  respect 
to  the  air 

ft/sec 

VC  (3) 

Vci 

Relative  velocity  of  the  c.p. 
of  the  "i"th  link  with  re¬ 
spect  to  the  air  ;(i-l,2, 3) 

ft/sec 

VG 

V 

g 

Gust  velocity 

ft/sec 

VGH 

VgH 

Component  of  gust  in  Y 
direction 

ft/ sec 

VGS 

Vgs 

Component  of  gust  in 
direction 

ft/sec 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

VGV 

Vgv 

Component  of  gust  in  Z 
direction  B 

ft/sec 

VHR ( 4 ) 

VHri 

Relative  wind  velocity  of 
"i"th  link  hinge  in  the 

Yg  direction , (i=l ,2,3,4)  . 

The  first  hinge  is  at  the 
winch. 

ft/sec 

VS 

Speed  of  sound  (not  used 
in  proaram) 

ft/sec 

WG  (8) 

An  array  of  eight  variables 
signifying  gust  velocity 
acting  on  balloon  and 
corresponding  to  TTG  (8 ) 

ft/sec 

WW  (8) 

An  array  of  eight  variables 
signifying  steady  state 
wind  profile  and  correspond¬ 
ing  to  AALT ( 8 ) 

ft/sec 

vw 

V 

w 

Steady  state  wind  velocity 
at  ZB 

ft/sec 

VWC (3) 

Vwci 

Steady  state  wind  velocity 
at  Zc.  ,<i  =  1,2,3) 

ft/sec 

VWH ( 4 ) 

VwHi 

Steady  state  wind  velocity 

at  ZcHi  '{l  =  lf 2'3'4) 

ft/sec 

WB 

WB 

Total  weight  of  balloon  - 
includes  balloon  and  bridle 
material,  and  enclosed  gases 

lbs 

WPL 

Weight  of  payload  located 
at  bridle  tether  connection 

lb 

WTC 

Weight  of  tether 

lb 

XB 

XB 

Lateral  displacement  of 
balloon's  c.g.  in  the  inertia 
reference  frame 

ft 

XBD 

XB 

dXB/dt 

ft/sec 

XBDD 

XB 

2  2 
d  Xg/ dc 

.  2 
ft/sec 

XBDDR 

Lateral  acceleration  of 
balloon's  c.g.  relative 
to  air 

,  2 
ft/sec 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

XBDR 

Lateral  velocity  of  balloon's 
c.g.  relative  to  air 

ft/sec 

XCD ( 3 ) 

X  . 

Cl 

Lateral  velocity  of  c.p.  of 
"i"the  link  , (i  =  1,2,3) 

ft/sec 

XCDD ( 3 ) 

r* 

X  ci 

Lateral  acceleration  of  c.g. 
of  "i"th  link^(i  =  1,2,3) 

ft/sec2 

XCHD ( 4 ) 

XcHi 

Lateral  velocity  of  "i"th  link 
hinge/ (i  =  1,2, 3,4) . 

The  first  hinge  is  at  the 
winch 

ft/sec 

XI  (3) 

Pitch  angle  of  "i"th  link; 

(i  =  1,2,3)  (positive-rotated 
up  from  the  horizon, 
clockwise) 

rac 

XID (3) 

• 

«i 

dq/dt  ,(i  =  1,2,3) 

rad/.'ec 

XIDD (3) 

• « 

d^./dt2  ((i  =  1,2,3) 

rad/sec2 

XIDDDE (3) 

h 

d2^./dt2  ,(i  =  1,2,3) 

deg/sec2 

XIDDEG ( 3) 

• 

h 

d^/dt  f  (i  =  1,2,3) 

deg/sec 

XIDEG (3) 

h 

Pitch  angle  of  "i"th  link, 

(i  =  1,2,3)  (positive-rotated 
up  from  the  horizon, 
clockwise) 

deg 

XID2 (3) 

• 

52i 

(i  =  1,2,3) 

rad2/sec 

XIODEG (3) 

^io 

Equilibrium  pit^h  angle  of 
!,i"th  link;  (i  =  1,2,3) 

deg 

XPLDD 

*PL 

Lateral  acceleration  of 
payload 

,  2 
ft/sec 

YB 

YB 

Down  range  displacement  of 
balloon's  c.g.  in  the  inertia 
reference  frame 

ft 

YBD 

yB 

dYg/dt 

ft/sec 

YBDD 

<  * 

yb 

d2Vdt‘ 

ft/sec2 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

YBDDR 

Down  range  acceleration  of 
balloon's  c.g.  relative  to 
air 

ft/sec2 

YBDR 

Down  range  velocity  of 
balloon's  c.g.  relative 
to  air 

ft/sec 

YCD(3) 

Y  - 

Cl 

Down  range  velocity  of  C.P. 
of  "i"th  link  ,(i  =  1,2,3) 

ft/sec 

YCDD (3) 

9* 

Yci 

Down  range  acceleration  of 
c.g.  of  "i"th  link. 

(i  =  1,2,3) 

ft/sec 

YCHD ( 4 ) 

YcHi 

Down  range  velocity  of  "i"th 
link  hinge  ,  (i  =  1,2, 3, 4) 

The  first  hinge  is  at  the 
winch 

ft/sec 

YPLDD 

•  • 

ypl 

Down  range  acceleration  of 
payload 

ft/sec2 

ZB 

ZB 

Altitude  of  balloons  c.g.  in 
inertia  reference  frame 

ft 

ZBD 

• 

ZB 

dZB/dt 

ft/sec 

ZBDD 

•  • 

ZB 

2  2 
d  ^B/ du 

ft/sec2 

ZBDDR 

Vertical  acceleration  of 
balloon's  c.g.  relative  to 
air 

,,  ,  2 
ft/ sec 

ZBDR 

Vertical  velocity  of  balloon's 
c.g.  relative  to  air 

ft/sec 

ZC  ( 3) 

Zci 

Altitude  of  C.P.  of  "i"th 
link  ,(i  =  1,2,3) 

ft 

ZCD ( 3) 

• 

Z  . 

Cl 

Vertical  velocity  of  C.P.  of 
"i11  th  link  ,(i  =  1,2,3) 

ft/sec 

ZCDD ( 3) 

0  * 

Z  . 

Cl 

Vertical  accleration  of  c.g. 
of  "i"th  link  ,(i  =  1,2,3) 

ft/sec 

ZCH ( 4 ) 

Z  ... 

Clll 

Altitude  of  "i"th  link 
hinge  t (i  =  1,2, 3, 4) 

The  first  hinge  is  at  the 
winch 

ft 

A  - 1 8 
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STANDARD  DESCRIPTION 

Z  Vertical  velocity  of  "i"th 

link  hinge /(i  =  1,2, 3, 4) 
The  first  hinge  is  at  the 
winch 


UNITS 


ft/sec 


•  « 


Z 


PL 


Vertical  acceleration  of 
payload 


ft/sec 


2 
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SECTION  II 


LONGITUDINAL  EQUATIONS  OF  MOTION  FOR  SMALL  ANGULAR  VELOCITIES 


First  consider  the  longitudinal  equation  of  motion.  The  dependent 
variables  (see  Figure  1)  are  0  (pitch  of  the  balloon)  ana 
(pitch  of  the  "r"  the  link) ,  where  r  is  a  particular  link. 

The  appropriate  equations  of  motion  in  the  Second  Scientific 
Report,  Appendix  A,  Section  II  are  (80)  and  (97).  Since  these 
equations  were  derived  in  three  dimensions,  thay  contain  the 
lateral  degrees-of-f reedom,  ip  (yaw  of  balloon),  <J>  (roll  of  balloon), 
and  ar  (yaw  of  "r"  th  link).  It  will  be  assumed  that  the 
balloon  remains  in  equilibrium  in  the  lateral  degrees-of-freedom. 
Therefore,  \p  =  <t>  -  ar  -  0.  With  this  assumption  Equation 
(80)  and  (97)  are  rewritten  as  follows: 

f.  ce  rc-tiff]}  &  y 

-  Fe  (i) 

i  mj’if'i clua  A'L 

t  dt-t)  fiJf  ce  +c£$e\[s&ce*c£se] 

**4*/  t~i 

^])  1  h'gse  tt£c&\[-s£  se+cf^aPi 

-Rf*jte[-stse+ctce>~]y <2) 

where  S0  =  sin  0  and  C0  =  ccs0  ,  etc. 


A-2.0 
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FIGURE  I  -  BALLOON  TETHER  MODEL  IN  LONGITUDINAL  PLANE 

In  order  to  solve  Equations  (1)  and  (2),  it  will  be  assumed 
that  N  =  3.  Rewrite  Equations  (1)  and  (2)  with  this  assumption 
and  substitute  in  the  trigometric  identities  for  the  sums 
of  angles. 

AMJfaU$+e)-4l¥Jlti„c(fl+6>'i)t 

*  • 

ilMi'Jfastfse)  +*)1  6  ( 3 ) 

+L  yj/tJt„$(fi+6>)-4i,//ij*  dt+e)l  fj  =Fe 


For  r  =  1, 


Ui  Mi,  5 

*f\r?  t/WFt  tHiS'Cf,-*)  MvL'lfaefiZ 

’  0 

*0 

Al^nUH J4sfe) 

For  r  =  2, 

[,iU sY/a  fsS e 

♦♦ 

,/[(!?  iMidtseM*)  ^MrddfsSS, 


»  • 

,/[  f  w*,,  +/>i  s*HL*)  ^  C2( 


For  r  =  3_ 


Uifo.  sifsd-MJt,,  c(.1t>*)\e> 


A-d.1 


l’l(  x*<*A)  cLfi-fi)  *A$ifi*0)$(fi+0)*Mvc(fi+6)  fj 

7  +»•„  t  A  5  Y£  *6)  t  re)}  'i 

Sj  (6) 

The  generalized  forces  (F^ ,  Fy ,  Fp  ,  )  were  derived  in 

Section  III,  Appendix  A  of  the  Second  Scientific  Report*. 

The  following  necessary  equations  are  taken  from  that 
section. 


^  ^  $  RaS  R-iiy  T  ^ty/t  “  fiHA  t?)fi  -1 

+  \~LsRjt  ~^/f  ^  ^GB  ^ 

f1,  'lU-W,  'hvt-/wn%  -bnv  £/, 

+!-/;„  -i*j JLst'Jc,  [  >  £>«,<: JJS  ] 

'jtllDctjS2?*  Dchi  $  fj 

't>^c‘?,)cl£-f,)l  (8) 


*Reference  2 
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(9) 


'[ L  rWt  +Ftu/I  ^-bfLV  cf»  +i~FtM  *AuJ^s Fi 

* 

-Alfe„S ^i-DcnC’Zl+UcS'f,  ^D^C1fJ]c(f>-f,)} 

Fft  '[  U '  K'i  f  F»„  ‘j  -Dnii  ■  T?  f\^  C^3 


♦  irFu*  -/>*  J/s  'FrJ.  £&■  * 2-£  ^ ;>6 ] 

(10) 

/>* '  DR/Urt  ( ’cf9)  +  L  IF  T8  (S  fa) 

(11) 

Fm  -  LIFTS  (C  h)  '  DR  At  t($  Fa) 

(12) 

i >to*S‘f.s,[£i,  *  ciie(^)F\ 

(13) 

line  *j,s,l£ »  +  ^»(^)<?] 

(14) 

/We.  -  y*  1  t  ^*hi  J 

(15) 

Dc^  ;  i  5tA-  £dc  Vc<  {  V^ci  ■*  Vu  ) 

(16) 

Dct/y  -  *3  /?4  Sc*  ^OC  14<  (  &*) 

(17) 

The  drag  on  the  payload  is.  neglected. 


bpt.H  '  tipLi, 
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(19) 


C  -  J -/  \ 


h  --  i  4  14,’ 

(20) 

Van  -  (yb  +  Vw  +  Vfn)  + 

(21) 

Xg  '  \m(iB’“Vftf)/(Yg  + 

(22) 

e-h 

(23) 

See  Figure  1 

In  order  to  show  the  load  distribution  along  the  tether 
at  any  time,  the  tension  in  the  tether  is  calculated  at  the 
end  of  each  lirk.  This  is  done  by  summing  forces  in  the 
vertical  and  horizontal  directions  at  each  hinge.  Before 
siamming  forces,  consider  the  inertia  forces  acting  on  the 
balloon.  Figure  2  shows  the  balloon  and  with  the  apparent 
inertia  forces  acting  along  and  normal  to  the  centerline. 


MA2gce-yebe) 


FIGURE  2  INERTIA  FORCES  ACTING  ON  BALLOON  IN  LONGITUDINAL  PLANE 
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(24, 


^Fi My  (iB  CG-vBu>)t&  +MJ'/,  ce+Z$Se?)se 

=  Hflt/^VrAsV]  +  sete 

\  Fy  =  M,.  liB  c&  +  2b  s&)  c&  -Mv(ia  C&  -9#  S&J  se 

~  ^a  LA)i'>My'}s8^8  +  y^I MyS  &tA1lC  s\ 

The  summation  of  forces  acting  on  the  balloon  are: 


ih  '  Ls  +  f bva  ~Wb~  ~Fbv  (26) 

^.Fy  -  Ink  ~~  f~an/l  (2/) 


where:  T_.t  and  T  are  the  vertical  and  horizontal 

bv  bn 

tension  forces  at  the  bridle-tether  confluence  point  respectively. 

Note  that  the  payload  forces  are  not  included  in  these  summations, 
but  will  be  included  shortly.  Equating  Equations  (24)  and 
(26),  and  Equations  (25)  and  (27)  gives  the  following 
expressions  for  Tgv  and  Tgg. 

t$¥'  ls  +fsva  -wB  - is  l rAh s V] “ yB  se c&  (28) 

Jbh:  Fbba  *  inlA'Alylsece  ^b  I MyS'e  +A  cz&]  (29) 


Equations  (28)  and  (29)  give  the  tensions  on  the  tether  at  a 

point  just  above  the  payload.  Figure  3  shows  the  payload 

with  all  forces  acting  on  it. 
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FIGURE  3  APPLIED  FORCES  ACTING  ON  PATLOAD  IN  LONGITUDINAL  PLANE 


Summing  forces  in  the  vertical  and  horizontal  directions  gives: 


*  », 

'™n  *  ei  '  Tsy  -/WfiL  j  -Fc 


VN 


(30) 


^  F  '  'if  i  ~  Fchn 


(31) 


where  F cv  and  FCH  are  the  vertical  and  horizontal 
components  of  tension  at  the  top  of  the  tether  below  the 
payload.  Substituting  Equation  (28)  into  (30)  and 
(29)  into  (31)  and  solving  for  the  vertical  and  horizontal 
components  of  tension  at  the  top  of  the  tether  results 
in  the  following  equations. 
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Fcvh i '  Is  Ft?#  W's  —  f  V  tMi,  S  Vi 

-y«  [/#.-/«,.]  C(p  ~~/yy] pL  -  />^? pl  ~ipt. 

fcn*z  Fma  +  ?#  \.rfh'Mv\  &+-/V1lC  &~1  Vpl 

The  totai  tension  at  the  top  of  the  tether  is: 


(32) 

(33) 

(34) 


Figure  4  shows  the  applied  forces  acting  on  the  "N"th  link 


FIGURE  4  APPLIED  FORCES  ACTING  ON  "N"TH  LINK  IN  LONGITUDINAL  PLANE 
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The  aerodynamic  forces  Dqjn  and  °CVN  are  the  horizontal  and 
vertical  components  respectively  of  the  normal  aerodynamic 
force  (Fcn,^  )  acting  at  the  c.p.  of  the  "N"th  link  due  to 
the  relative  motion  of  the  nN"th  link  with  respect  to  the  air. 
This  relative  motion  takes  into  account  both  the  steady  state 
wind  and  the  horizontal  and  vertical  inertial  velocities  of 
the  link.  Expressions  for  Dchn  and  are  given  in 

Equations  (16)  and  (17).  The  normal  force  acting  on  the  link 
is  given  by: 


~  Dc/t„  $  fu  +  Det/»  £ 


(35) 


Summing  vertical  and  horizontal  forces  on  the  "N"  link: 

^  Fi  "  /M  2.V  “  f~£VN  ~  Fcn^  C  f/j  ~  —  Fc  V  M-!  ( 36 ) 

-  /tV  Ycy  *  ~  Fth¥  ~  Fcti"  $  f«<  +  Flhi  ^ 

Solving  for  the  tension  forces  at  the  bottom  of  the  "N"th 
link  gives  the  following  equations. 


Fcvtbi 

■  Few 

Cfj¥ 

(38) 

Fchm-i 

Fen 

tfi cNh  S&  *  /fr  YcN 

(39) 

F :uA-t 

--  fc* 

-  J 

(40) 

Fc  Ha-/ 

-+Fcna  SfA  t/w  Yca 

(41) 

x, 

(T 

»  I- 

1  f  1 

./  *  Icha., 

(42) 

The  tension  at  the  winch  is  found  by  setting  v  =  1. 
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SECTION  III 


NON  LINEAR  LONGITUDINAL  EQUATIONS  OF  MOTION 


Section  II  of  this  Appendix  contains  the  longitudinal 
equations  of  motion  of  a  tethered  balloon  if  angular  velocities 
are  small  and  products  of  angular  velocities  are  negligible. 
This  assumption  was  made  in  the  derivation  in  Appendix  A  of 
the  Second  Scientific  Report*.  It  is  now  desirable  to  examine 
the  equations  of  motion  when  large  angular  velocities  are 
allowed.  As  in  Section  II,  appropriate  equations  will  be 
extracted  from  Appendix  A  of  the  Second  Scientific  Report*. 
These  equations  are  (41),  (79)  and  (96).  Equation  (41)  is 
the  total  kinetic  energy  of  the  system:  Equation  (79)  and 
(96)  are  terms  that  are  needed  for  Lagrange's  equations. 
However,  these  equations  are  functions  of  both  the  longi¬ 
tudinal  and  lateral  degrees-of-f reedom.  Therefore,  before 
rewriting  Equations  (41),  (79)  and  (96)  of  Appendix  A* 
set  all  lateral  degrees-of-freedom  equal  to  zero 

(^  -  <t>  s  ar  s  0) . 

The  total  kinetic  energy  of  the  system  in  the  longitudinal 
plane  is: 


/a  -  f 

HA ( V*J>)  HA i*  ** 


(43) 


where 


(44) 


*Reference  2 
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(45) 


r-n  i  % 

Vs- J$  £  Z 

Now  write  Equation  (79)  of  Appendix  A* 

*My{/2  i  c  (&+ 1 ) '  e  U- it)*, }  +Ijs& 


Rewriting  Equation  (96)  of  Appendix  A*: 


ff  =  |  <* r\ %  i  c(t-t )  *ilc(r.  -fj 

fMt-t)  *iii  ,MnJV  l  C  (£-t  )} 

a:  ‘  « */ 

+  Ml \j'f.  1  $  ($*£)  t) 

A  i  * 

+A)Jj2^  jj.  Ue-tfx)  \  c(&+f«) 


Langranye’s  equations  for  the  system  are: 


'  fr "  F* 


sL  (Hr  )  -  2l  Ft 


*Ref erence 
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First  derive  the  terms  in  Equation  (49).  Using  Equation  (47) 
the  following  term  is  formed. 


i(6'6)d#£)\*b.A  &- 

*  -i  *  <-/ 

+  1«  l 9  (51) 


From  Equation  (43) 


From  Equation  (46) 


7%  f,  <'if* 


A  -ii 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 


w “  r  ,  . 

£  -4  =  4.(9 

(58) 

V*4?'*jdr  £de*i)-L^e 

cO  *  - ' 

(59) 

(60) 

06 

Substitute  Equations  (53)  to  (60)  into  (52): 


de 


i  i  b»t)  eK  4.  6*AicUrfJ- e] 

*H  ’  •<-■/  • 

«-#  r  /:/ 


(61) 


4?  -  u-a]  r/£  i  utht) ] 

v*-' 

J  “  A, 


(62) 


Substituting  Equations  (51)  and  (62)  into  Eqation  (49) 
yields  the  final  nonlinear  balloon  pitching  equation. 
N^te,  Fn  is  the  same  for  both  the  nonlinear  equation  that 

U 

follows,  and  Equation  (1)  which  assumed  small  angular 
velocities.  The  expression  for  F0  is  given  by  Equation 
(7). 

» Jh.  1 1  s te*£)  -MJL„ i  i  dU\  *k«L  <An‘  <L.k 

a -I  *  -i:i 


-  1  i  stet&ii?  £  c(0+&) ] 

A-! 


\-t  4-i 


h 


(63) 
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Now  consider  the  tether  equation  of  motion  (50).  From 
Equation  (48)  it  follows: 

_  M  /*-v 

i  cit-t)  -ltd  -i)  *  it- 1  j 

v*  J  -n:~n  >i/  V-l 

4  Idt-iHid-i)  j  (t  fi  w/i  fit-  c «- 1 ) 

i(UMi-i)  t  di-Q 

a  ;  i  *i~t 

-i  id'tuit-t)] 

A 

\l£‘t  i  s  te  *£•)  d  '1 1  Iti.  t)  de>d)  dku.  e\dt»£) 

+lft  Ui^DdK^fs)  U.iAc(e»D\ 

i-l 

>,d\d'Vide.t\-l'l  iUd)sdd)d^  file  le»L) 
-  lj2 1  £  *  j[<?  + 1]  s  fef/J 


(64) 


From  Equation  (43)' 


fclM  v*-  If" )  '***  ( v*  j£/L] 

’i»)  +A  (VU*-&g) '  Mt ) 


(65) 
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From  Equation  (44): 


1  -If*  £ 

A.  'y/l 

H  ' 

i 

J  -yi  £ 

z«f/l 

(66) 

1 

0 

/*i  </t 

(67) 


From  Equation  (45)' 


iit-n 

it 

From  Equation  (4^./; 


(68) 

(69) 


'  L  '  -J  Ll-  c  fAca  t  i  &  $<?] 
il*t  *-jtLlctse>sf«ce] 

dl 


(70) 

(71) 

(72) 
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Substituting  Equations  (67),  (69),  (71),  (72),  (54),  and  (55) 
into  Equation  (65)  yields. 


& -  jm  {f  [  f  i  i  ut- -t)  1 51 a  -£i\ 

-if  ii s(t-t)}  -^nrt iis(t-t) 


if-i 


-All  c ie* t)Ul  t  ile*£)  -  IL~  el 
Lsle+t. Jfi7^  ic(e+  %)-/?)*«  ei] 


Substitute  Equations  (64)  and  (<  ■)  into  Equation  (50)  to 
give  the  final  nonlinear  tether  pitching  equation. 
has  already  been  found  in  Equations  (8),  (9)  and  (10). 


fmjt\ r  l  dll)  l  dlfj  wj/i|  i  dl-l)  4 1  j 
-m >,/!  idl£)-/KUi 

in  <si 

-Mv\ri  l  c(s*l)  -J  (j.„  tt\c(e*1„) 


=  t**s{5  ils(i-t)-iils(L f.)l  w'f  i-sti-t) 

--A+I  <x/ 

+/MPLJr  '{ -Ale'[i'Us<i)d'»t) *[**£,) 


A, 


(73) 


(74) 


i-i  (e+i)sfo+£)  c(e+l) 

A  l  eue+t)5(e  ^)  -JLi f,-*.  el$(e+fj  1  + 

*si 

Now  assume  the  tether  is  composed  of  three  links  (N  =  3). 
Rewriting  the  balloon  pitching  Eq ration  (63)  gives: 

{ he  iu  +/t1*  $1*  ?  & 

•  » 

+1[m,  4*  5  C  (*+  fJl  /» 

3  i 

islsutn  icte+f,)} 

A-l  •* 

t  cie*t)-A\£  ft, Jz  i'slp+l)  *Fe 

izi  *  =< 

For  N  =  3  and  r  =  1,  Equation  (74)  becomes: 

UJKa^  >(0+f,)  '-MvJtlf*  C  le+f,)]  9 

+iVf  *'"»n  *Mk$x(e+t)  +M¥cx(e+t)}  f, 

tf-fj  s(e*f>is(6»f,)  iM¥cl#+&)  cM,)}  ft 
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»/ ’{[ ?*■««] c(frf,J'MiS : letf,) *A1,C (e*f,) c(etf,)}  f, 

3  J 

+l/Mt,'/rtJjis(e+f,)l  iec(e-Z)+c(&+ft)l  % &s(&+t)} 

+M¥£2c(&+f,)?  i*s(e+Z)-M^/ls(6+f,)%  f/c(&+&) 

+  -!  <=l 

+Fft  {76) 

For  N  =  3  and  r  =  2,  Equation  (74)  becomes: 

\AJUa*  s(e*f>)-M¥J  ftj„  c  Le  +  fy)]e 

+f{  l  c(ft )  +Ks  (&  *■  f,)s(e> + t)  +M¥  c(e+f)c  (9  *■  fj  J  f 

t /*{  t  sUe  +f*)  ?* 

t f\  lr*l+/*v)  c  ( f»-  f3)  t  AS  (&+  fi)s(&+  fs)  +Myc(p+f,)c(e*  fjJ  j} 

=^Jl(  i  t  \  ( frt ,)  4  f'sfofj  Wfifl  I**  ( t  ~  S9) 

*Z1 

AA^JJ2{s(eftf9)i  iecfatl+cfatfjji  %6*(e+t)} 

A*  I  4  -  i 

a  r/  A  r/ 

*-F1f  (77) 
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For  N  =  3  and  r  =  3,  Equation  (74)  becomes: 


c  (<?+&))  e 

t/f l^+Mfhdfrf))  +AUe+$)sCe+1i)  +MvC(ai-fl)cie>ffJ]  f, 

y*  f>) 5  feifMMrCU &) c L 9+ fi)]  & 
t/*{f  ***h  +AS1(e+fi)+/rtl,c*(<?+fi)]fi 


+  lA~/M‘-U1U(e+fj)5  ieUfi+()+c(e+f,)Z  &»s(*+£$ 
+A*J\ c  le+fjt  i*s(p*t)~#i.jta3te+fi)f  i’t  (e+&) 
-JlA  &*.  cU+%)  +/tfr/(j+,s(s>+fij}/?*  +  Ff} 


(78) 
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SECTION  IV 


LATERAL  EQUATIONS  OF  MOTION  FOR  SMALL  ANGULAR  VELOCITIES 


Sections  II  and  m  considered  the  equations  of  motion  of  a 
tethered  balloon  in  the  longitudinal  plane.  This  section 
will  consider  the  equations  of  motion  of  a  tethered  balloon 
when  the  longitudinal  degrees-of-freedom  are  fixed  at  their 
equilibrium  value  and  the  system  is  allowed  to  move  only  in 
the  lateral  degrees-of-freedom.  Figure  5  displays  a  front 
and  top  view  of  the  tethered  balloon.  All  angles  are  shown 
positive.  The  lateral  degrees-of-freedom  are:  ip  (yaw  of  balloon), 
4>  (roll  of  balloon),  and  ar  {yaw  of  "r"th  link). 


. '  KJ05C  OF 


top  virw 


FRONT  VIEW 


FIGURE  5  -  BALLOON  TETHER  MODEL  IN  LATERAL  PLANE 
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The  appropriate  equations  of  motion  are  Equations  (73), 

(86),  and  (108)  of  Appendix  A  of  the  Second  Scientific 
Report*.  These  three  equations  will  be  rewritten  below 
remembering  that  0  (pitch  of  balloon)  and  z,  (pitch  of  "r"t 
link)  are  constants.  The  first  equation  is  the  balloon  yv:v 
equation  of  motion. 

N 

Mi  {T£'i[w  ktW  s*4Ltc'r)  +  S%ct(s+s0C't-  c4  5  •rj-Saj  5 £  (-Sjcfj  j 


»■[  Rj*  (-vcec*)  -4*  (-fsj+t)] flu.  s^S 

4rtM%$Uz(cosv)+S£6t(cecd-s(Gst(seMlb'(*s*&$1hrtc*l 

4  SI 

+41 M  £  V}  [ c<£  (~  Sfi  ty) + S£  c  i l~  c4  sf  c<r  -54  w)  -  $<£  5  f  ( £  ,*  £sj\ 

aU 

+  1  'tjm  ( V  S4  £#)]}  {-#),„  $4c  9} +4 1-  lye  M  *  l  v*  e  CG  C4  ] 

f  I  L'dllx,  S34  +  liiC  If)  tSX&lyg  ~$2&  Ctplyit]  -  ’f  (7 


The  second  equation  is  the  balloon  roll  equation. 

Mi Ic(%(s4>s0$</+c#c<i>)+;cc&(s4> $ec*-c4>$<f)-${Z5f  (-S4ce 3 

+ 1  £,*,(-&  ct)  -  4*.  (-  ss?  +  $  )j  ]  I*  4*  J 

+  </>  l  lyt]  +  ^  l~Iu  $6  r  I/*i  C&CQ ]  -  4  ( 
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The  equation  for  the  tether  link 


given  as  follows: 


yaw  equation  of  motion  is 


i* 

A-l 

N 

j' 

i  crAcrn  *s c  *6  c  (i-f,\\ 

Ucfasesv***  Cf)+3G  ctlstsGc*  ~aw)-S(TA  5fA[^<pce)} 
('^cecf)-^^  C-v$e+$jl]'jjlctz(s4>s&sip 


*  &  c-t  l$?  5o  opsyj-w  ;>£  i^c&j]} 

[CZ  (ces</)  +s%  c £ l c& cv) ~ $or  <,£.[> e) J 

]{JlcGlcestr)*SGi  cZ(cectf)-SK  $£  ($*;]] 
rMvU$ti Uk  r$fc *)*■$(% c£(-c<t> sec*  -sfst) 

~s<?32  teteft+i'Rfm 

*  s&  c&  I't4>$ec#-s4>str) $  £  (c^c»j\)-  F 
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In  Equations  (79),  (80),  and  ( 81 ) t  there  are  several  trig- 
nometric  combinations  of  the  balloon  angles  4>,  0,  d>.  These 
terms  are  the  elements  of  a  matrix  (Equation  (51))  given  in 
Appendix  A  of  the  Second  Scientific  Report*.  The  matrix 
(^j)  gives  the  relation  between  the  balloon  unit  vectors 
and  the  inertia  frame. 


V 


5  4>  se  Stf  +c4  c* 

s4se?ctf'C4s* 

-s4ce 

ce  $<r 

ce  cy 

$e 

-cfsectf'-s^sv 

c4  ce 

The  balloon  yaw  Equation  (79)  is  rewritten  as: 


i'l 

+/H \jl l ctr  ilt  t  bcr;  c £  $£  6j][- 

«•*/ 

r  N 

Ml, u £  ojl  car-  cit  +  c¥,ei3-s<r:  s &  *',3]  f  fy.*,  c,*} 

A  ~l 

£ji  ^li]  f  ^  *»  f  Pj*,  £‘i\Mv 

+  ce[h*$*4>  *ha  CV]  +Sx&ly9  -S20C4>ly£6  f 
*  s  *  lv38  CfiCip]  ' 


(82) 


(83) 


The  balloon  roll  Equation  (80)  is  rewritten  as: 
*Reference  2 
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Ai/f 

<» 

I~/;*  £jj~*  &,*,  ^w]['/ifc/i,l/^i  -  Iy&  S&  +ly*e  C^4>] 

+1-Yg$ 


(34) 


The  tether  yaw  Equation  (81)  is  rewritten  as: 


+s<z  s%  c{fn-t)]+j[<L[cfc 
+  tJ*T  %lcG£G+SGSiGc(t-Z}j]+¥Jti 

4  »/ 


CtfT  -*S<E,  50TC 


(£-&)]} 


?r 


ia£  £  *507S0:^fA*£)] 

^  *) 

i  L&T  £„  ♦54F't£  £'*  ~$tr:5t£' ill*K  6»  +$<Z£&*u‘$KSt  £<i\) 

<*•/ 

*Mu{jlf&[cQl„  i-saatfiCu  -Saji  £  £»,]U£^+-C£c&  t*  -*<£  sl^«8 

A~l 

+M,\f'LvAw£j>+s<r;L i  eis  'ior ij  ^ i*sizcC £»-<>& st<£»l} 

A :> 

It/fk  [-#i„  iii  +Ru  £«McK  t„  i^cIa  -so;  St  6,)] 

A  l~$tsm  £ji»  £*j  *'£  £jjl  tn  ^/jt] ^3/ 
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In  order  to  solve  Equations  (83),  (84),  and  (85),  assume 
^  ~  3.  The  balloon  yaw  equation  is: 

£{Mi  [c^i,  *  W  Cf,  6n  -SGif,6i)]l-ltjtm6„+fiAm  eti\ 

caifr  nsicf,  urscr,  sf,  £„]  [-  £u.  e  M] 
+Mt,lcreti,+st;cf,e}i-srtsft  £*]$c> ; 

d{toX cn  6»  +$o»  cfxtn~  Sft 

*M,lcrt ta,  *  s <s  c£  c»  - s tj s £ £ ,3~] 

^  f  6*  rMcf,  £ir  M  S  f,  £  J3]  [  •  w  £,$  71 

'JlMilcq  £.  1 50}  cf3t„  -  5rM3  £,}]hj.„  £„  f  £Zi] 

+MAc<ri  £„  t5oj  cfj  6z s<r;  s£  £„][-  &,»  £„] 

It<ruit-H><r,cf3  £u  -  sjj  s  f, £„][  /?  ?>  £ ^  j  ^ 

I  £jj  i  ft  A  A,  £aj]  +-M).  £'yjl  +4^1 
^Vl-Xn  t  +SVlyp-  Sl&  c4>  ly#a^  V 

+^A*>  ^»}ll ~fo*i] ~lyg$&  £6  teipS  41  ~  F^> 
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The  balloon  roll  equation  is  ' 


j\AUicr,6„  +s <r;cft  etl  <r, 

ij  \M%  l  cr3  L„  i  s<ra  cfz£„ -Wish  £ii)l-Kj,~]} 
tjfyslco]  €„  +S0]  cf^n-%0}  s/j  <r3 

AtisVtj*.  -lye  W  ^lnHC,Bc^ 

+iy e}4>  -  F* 

From  Equation  (85),  three  equations  will  be  generated 
(one  for  each  link)  by  assuming  r=l,  r=2,  r=3 
consecutively. 

For  r  =  ly 

Jl  j  ~J~  +  /W Pi  i-Ais  L  C0~/  (n  +  S^T  ft  £  a  ~  ^  1 

+^i  L  cr,  it,  t  S07  cf,£n  -W  s  1 l+M¥  [  cr,  £„  +  S  ^  ft  £  H  ~  S  ft  6  J]  C 

tj*{  i  ur%  ctr,  a  S<7i  Scg  c(  ft 'ftj]  [  TT*  +/nfpi\ 

+Ah[c<r}  £„+$<%  eft  £u  -$<r,  sfre^U^e,,  ft  (=n '  5*7  s  ft  Oil 
^41  [  1 03  £a/  +  S(Ta  c  ft  £  a*  -$0S  S  ft  6«l[c*70/  rS^C  ft  £„-  S(T,  5  ft  L zjl 


(87) 
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■  IW*9  JU.MNUU  ******** 


!L»JJ|p.iJ 1 


+/M  v[an  6),  r  5^  £&  £u  -  5flS 

1  ctr3  c<r,  +  ss}  $67  c(frfi J]  f  ^  *v»v) 
ttfilcni*  tSff]cf36u-'S<fjS'fl6»l  tc<r,6„+SJ;cf,£,j-$ff7$f,£,t] 
tfiUejtj,  +SG£f,£»-sc5sf36JtJ][c<r;6j,+SJ7cf,£ti-S<rj 
+A„Urs£Jt+$o 5  eg  fja-SflSSgtjilU^ii,  'WeftS*- SQ->s1,£-n]S% 
t 'tfj/tm  £»*S<n  cf,  £/2~$0/Sf, 

+  t-tr,  £21  *“5 <77  C  f,  in  ~S(77  sfi 

■MWify*,  £j[cr,  £i/+S^£f,6i~S<r;5'g£jjj)  ^ 

£  6 1,  ~t$&7  £  f/  S&7  (88) 

For  r  =  2y 

f\Lar,art  +$rt<><rrL(t-f,)\l^r  **«»'] 

■*Milc(r,6„i$rt  cf,£n -  5*7  ytf  £,j]  [ Co*,  £ //  +5 ^  c f,  6*  - 5 ft  S  fx  £>t 1 
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t/  (  X*  f/r,;A  ^'/  £  /;  £>a"  5^,1  S&  £/jJ 

[iff  £»+!>«£&  ixt-S<n  Sf3  6i^\\mAc(Tx  6j,  +S<TxCf>  £s>S<rx 

tiX\'lCriC(Ti  +-50}S«  £"&,)][*?  +**/»/] 

*^s[<;q  ^,  +  SfiJ  c &  6,  -  s^5  $  fj  ^ilt  c*i 6"  * S<r3  C  f2  6,*  Sffl $&  £/j] 
+AI  Irn ix, ■* strj  r-fi  £xx  ~ $<r3  s  ft  6  A}]U^  **.  +*£  c£  6 ,  -  J  7i  sfx  £n ] 
+/lf,[c<>5  Cfa  iix'  S or,  SfUjj]  $  ^ 

£jj  £*3  Jf £« *  $0j  £fj  6/*  'S0»  S  ^/j 

tyflil'fo*.  £/jl  [coi  6i,  *$Ox  cfx£ix-S<rx$fl  £n~] 

MU  I  Rj~  £/*][c<rx  6h  t5<TaC  £  6ja  -SOiift  £si]j  £ 

t£  Ui  Vfij,*  Her,  e„  +$c>i  c  fx  i„-  $r. ;  $  A  £„]}  4  '  />t 


(89) 
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For  r  =  3/ 


/{lcr,cr~j  +$<r;$rj  c(frf, j]  [ ? 

+Ali[ar, i„  iSffscf,  £,t ~sr, sf,  6,^\ \ca j  £„  ^tr3cf}  &n  -$<r3 s  f3  £ti] 
iAUr,£3,  +s<r,cf, £n- S <n $  f, £«][<: r3  6i,*s<rjCf36»-s<r3  s f3  £*i) 
MyUfjtj,  ei2 ~sm  sf, <?dU<0  631+  so;  c  f3 ej,- 5^  s f3  r, 

*f  \ lc<ncrj  +Saz$<Tic(?r  fJJ 1 

+AU I cr3  64  *  sracft eu-  sai  sf,  £,3] [  c*i  t„  <-  stri c  f3£,3-Sg}S  f3 £/j! 

*AL  [&73£a,  tsos  cfUxa*  $<7i*  SfUdU*!  6u  +S05 ^  fi  txx-  50j  sfi  £ii] 
tM,[crxei,rf'r*cfa6i,-vi*T,et})[aii  £j»+*r9cfi£n- 5<*}S?3  £axl?  r2 
h/xi  T  +'*»»#«  t/^sf  c  7}  £„  +5(Ti  cfi  6,1  -  5c/}  s  fi  £> 3  ]l 
Uilcc 5  57}  cfi  £«  -5fi}5  &  &»3  M1,[«g  f  J.  *6<r}  c  /j  5  f4  6  J*}  r3 

t/'llji \~fjtH  633+^1^  6ii][<^}  £*  +scr3  cfj  6,2  S7i  Sf3  £73^ 

+/rill‘lh,£i3)Urt  6/  +str3  cf3  £xr  s*}  s%  £2J1 


A -49 


cfj  6to  -$tr3  sf3£*]]f 

*MAkl-&*\[cai  £„+$il  cf3£„-wi$  fa  f/j]}  4>  *  Frj  f90) 

The  generalized  forces  (F^,  F^,  FG  j ,  F^,  Fa3)  were 
derived  in  Section  V,  Appendix  A  of  the  Second  Scientific 
Report*.  The  following  equations  are  taken  from  that  section. 

Er FtSA  {(-(?„  ifaA  6n)b  t  ifaA  £*)& .  *«)**]  (91) 

F*Tlu/s  fay  "L(  fai\  £,}  +  I*  Fat#  (?A*\  6/  +  ^#0  ^ 9  2  ^ 

/>,  --i/Ur  H'r  6-7*  *  I?)?]['±<r7  SfJ  +L FSS/f][£<r,] 

4  P‘A  U<  c\]  4Dc»]  [cVx  (critv;  +$c;  so;  c  (frfjft 
*1  faille1*}  (c<n  ur,  t!>0}  Wt(frfa)/J)  (93) 

Fji'jflu-Wt  -l^n  fallen] 

^/p<w]Ua  tSIiSfiicffj-fJlJ  (94J 

Fri*j{lls'Vrl'*k*2)f] f‘*5sfiW WlM  AhsJk  cln]  (95) 

The  aerodynamic  forces  and  moments  are: 

^ 5  u5*  i  ^ f  Ly*i%) v  *  £f*ia )  ^  (96) 
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(97) 

*  ****  (%)*] 

(93) 

bcu  ~  ~  2  fa  $a  (oc  Va  Xa 

(99) 

fr-kf.V* 

(100) 

Vm  '  (*« -tys)*  ♦(>'«■*•  MO*  +  (i.-l^v)1 

(101) 

r8  ^'(Vy„) 

(102) 

(103) 

(104) 

As  in  the  longitudinal  dynamics,  it  is  desirable  to  know 
the  tension  in  the  tether.  This  is  done  by  summing  forces 
at  each  hinge  point.  First  consider  the  apparent  inertia 
forces  acting  on  the  balloon  along  and  normal  to  the  centerline. 

flu  -Mi  U'e  i„  t  y*  6u  f  He  (105) 

Ml  t  x*  U,  t  y#  £»  +  z*  £**1  (106) 

fj/  ♦  V*  £j3  +  2e  (107) 
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In  Equations  (105),  (106),  and  (107),  X*,  Y*  and  ZD  are 

W  x3  n 

inertia  accelerations  of  the  balloon:  A™,  and 

are  inertia  forces  in  the  positive  direction  of  the  balloon's 
lateral,  longitudinal,  and  vertical  axes.  are  elements 

of  the  matrix  defined  in  Equation  (82).  Summing  the  apparent 
inertia  forces  in  the  lateral,  longitudinal,  and  vertical 
directions  gives: 


s  £//dXe  * 


(108) 


4  Cxtdi'fB  ■*  £■  »s  t  (109) 

f  £*i#ya  +£u&»e  (110) 


The  forces  acting  on  the  balloon  above  the  payload  are: 

*  F$u~  T*s  (111) 

5.^-7 '$//"  Paha  (112) 

i*  Ls~\^t  * Fsv a  ~~Tbv  (113) 


Now  equate  Equations  (108),  (109),  and  (110)  with 

Equations  (111),  (112),  and  (113)  respectively,  and  the  following 

expressions  for  TBg  TBH>  and  Tgy  are  found. 


“Tg S'  F&tfF  [  £h  Cty$  *  £31 


(114) 
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(115) 


~f*V'  Ftvp^  L$"  “  (&1  Fix i  +6n&ys  +  CilClib) 


(116) 


If  forces  are  summed  at  the  payload  location,  the 
components  of  tension  at  the  top  of  the  tether  results. 


fix*  :  f*SA~  C  £h&B  +  £kt$yg  4  £)idiB  )  -/hlflXn  (117) 

Fen*  ~  FfiHA  4  ( in  $xa  f  €%%  l  £ja  i  ^pi  Ypi  (118 ) 

fcv*  '  Fwp  -f  ^  -  U/g  ~  l6n^Kt  t  i  L\}CU&  ^  ip I  (119) 

The  vertical  and  horizontal  forces  on  the  "N"th  link  are 
the  same  as  they  were  in  the  longitudinal  dynamics 
program.  The  force  in  the  lateral  direction  is  given  m 
Equation  (99).  In  general  for  the  "r"th  link. 

—  •  * 
ics.-,-,  -  h£^  -bc'M  -  s*n  Xco. 

FcHn-l  -  Fc*<t  ~t  Fcft^  5  £,  4,471  yc<, 

Fc\Jn-t  -  Fcvfl  ~  Fn„  c  tr.Wf  ~,y>1  2c* 

- - — - - — 

To-'  '  \(Fs /V/ )  (f £M-/) 

The  tension  at  the  winch  is  found  by  setting  r  —  1 


(120) 

(121) 

(122) 

(123) 
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APPENDIX  B 

DYNAMIC  RESPONSE  INPUT  AND  OUTPUT  DATA 
LONGITUDINAL  AND  LATERAL  DYNAMIC  SIMULATIONS 


CASES  1  THROUGH  40 


SECTION  I 
INTRODUCTION 


Appendix  B  consists  of  the  input  data  used  in  making  tee  computer  run 
cases  1  through  40, and  computer  constructed  plots  of  selected  output  variables 
for  these  cases. 

Table  B-l  lists  input  data  for  the  longitudinal  dynamic  cases  1  through 
23,  and  Table  B-II  lists  the  input  data  for  the  lateral  dynamic  cases  24  through 
40.  Tables  B-III  and  B-IV  list  the  balloon  longitudinal  and  lateral  static 
aerodynamic  coefficients  respectively. 

Computer  constructed  plots  (3  plots  per  case  are  included  on  each  page) 
for  cases  1  through  40  follow  Table  B-IV.  The  parameters  plotted  with  respect 
to  time  are  listed  below. 


Longitudinal  Cases  1  through  23 


Balloon  pitch  angle 

Link  3  pitch  angle  Plot 

Link  2  pitch  angle  1 

Link  1  pitch  angle 

Tension  at  bridle 

Tension  at  winch  Plot 

Balloon  altitude  2 

Balloon  range 

Wind  gust  velocity 

Pitching  velocity  of  balloon  Plot 

Wind  angle  of  attack  3 


Relative  velocity  of  balloon 
Lateral  Cases  24  through  40 


Balloon  roll  angle 

Balloon  yaw  angle  Plot 

Balloon  yaw  rate  1 

Balloon  sideslip  angle 

Balloon  lateral  displacement 
Link  3  yaw  angle  Plot 

Link  2  yaw  angle  2 

Link  1  yaw  angle 

Wind  guot  velocity 

Relative  velocity  of  Balloon  Plot 

Tension  at  bridle  3 

Tension  at  winch 
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Table  B-l.  (Continued) 


-56.6  3.7  -57.  5  ii.b  -5-4.7  8.7  -57.5  5.9761  49.0389  58.4903  70.9483 
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Table  B-II.  Input  Data  for  Lateral  Cases 


Table  ”~II.  (Continued) 


£ 
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.0433  11374.  14.7  164  4  154.4  1.94  0.0  -3.07  0.0  -.398  -.508 
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TABLE  B-IH.  LONGITUDINAL  STATIC  AERODYNAMIC  COEFFICIENTS 


TABLE  B-IV.  LATERAL  STATIC  AERODYNAMIC  COEFFICIENTS 


Sideslip  Angle 
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.  7260 
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-.  726 
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BJ  C 
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(Nominal)  n 
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193 
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1  0.  0 

.  193 

C 

y 

.  299 

0.  0 

-.299 
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.0014 

(Nominal)  n‘ 

C1 

002 
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.  002 
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502 
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.;  154 

(200%  Bottom  m 
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0.  0 

.  3065 

(300%  Bottom  m 

Tail)  C1 

0768 

0.  0 

,  0768 

C 

,  6b3 
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